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The energy and entanglement spectrum of fractionally filled interacting topological insulators 
exhibit a peculiar manifold of low energy states separated by a gap from a high energy set of 
spurious states. In the current manuscript, we show that in the case of fractionally filled Chern 
insulators, the topological information of the many-body state developing in the system resides in 
this low-energy manifold. We identify an emergent many-body translational symmetry which allows 
us to separate the states in quasi-degenerate center of mass momentum sectors. Within one center of 
mass sector, the states can be further classified as eigenstates of an emergent (in the thermodynamic 
limit) set of many-body relative translation operators. We analytically establish a mapping between 
the two-dimensional Brillouin zone for the Fractional Quantum Hall effect on the torus and the one 
for the fractional Chern insulator. We show that the counting of quasi-degenerate levels below the 
gap for the Fractional Chern Insulator should arise from a folding of the states in the Fractional 
Quantum Hall system at identical filling factor. We show how to count and separate the excitations 
of the Laughlin, Moore-Read and Read-Rezayi series in the Fractional Quantum Hall effect into 
two-dimensional Brillouin zone momentum sectors, and then how to map these into the momentum 
sectors of the Fractional Chern Insulator. We numerically check our results by showing the emergent 
symmetry at work for Laughlin, Moore-Read and Read-Rezayi states on the checkerboard model of 
a Chern insulator, thereby also showing, as a proof of principle, that non- Abelian Fractional Chern 
Insulators exist. 

INTRODUCTION 

Band theory formed the foundation of solid state physics for the past century. The theory analyzes the one- 
body motion of electrons through a crystal and obtains their wavefunctions as Bloch states dependent on a crystal 
momentum which takes values in a space equivalent to a d-dimensional torus (in d dimensions). Insulators hold a 
special place in band theory: at a first look, they are its most boring aspect, having a full gap between occupied and 
unoccupied bands, no low-energy excitations, and hence no interesting properties at zero temperature. This viewpoint 
permeated the physics of much of the last century, but is, as we now know, untrue. The space of the eigenvectors of 
the Bloch Hamiltonians can be thought of as a unitary matrix, which in the case of insulators exhibits several gauge 
symmetries related to permutations of the occupied and un-occupied bands, thereby creating a complex manifold. 
The Hamiltonians can be regarded as maps from the Brillouin zone (BZ) to these manifolds, which, as mathematics 
teaches us, can be nontrivial. Anytime an insulator with a nontrivial map (Hamiltonian) is placed in the vicinity of 
the vacuum, gapless edge or surface states appear on the boundary and cross the energy space between the valence and 
conduction band. One of the first examples of such behavior was Haldane's Chern insulator model on the graphene 
lattice 1]. This model exhibits the physics of the Integer Quantum Hall (IQH) effect but does not have an overall 
applied magnetic field, thereby preserving the translational symmetry of the initial lattice. The IQH does not require 
any discrete symmetry to exist. The field of topological insulators has further evolved to include other symmetries 
such as time-reversal charge conjugation and point-group symmetries. 

Most of the research on topological insulators has focused on one-body electron Bloch physics. While several 
interesting works on insulators with interactions exist [SHZ], most of them focus on obtaining a mean-field state with 
topological properties through interactions. Recently, it has been proposed and substantiated [5HTU] that fractionally 
filled bands of Chern insulators at 1/3 and 1/5 filing can exhibit an incompressible state in the same universality 
class as the Abelian 1/3 and 1/5 Fractional Quantum Hall (FQH) states. This state was dubbed the Fractional 
Chern Insulator (FCI). Its appearance is quite surprising, as the FQH and FCI problems vary considerably in several 
important ways, in the lack of a constant Berry curvature in the FCI, its large lattice filling, the absence of holomorphic 
and anti-holomorphic (z, z*) structures in the lowest band and the presence of Umklapp processes that favor density 
wave states. 

In a previous paper we have shown that the groundstate of the interacting 1 /3 fractionally filled band with Chern 
number 1 of the checkerboard Hamiltonian supports a 3-fold degenerate groundstate separated by a gap from the 
excited states and which is a featureless liquid with momentum orbital occupation number 1/3. We showed that the 
spectrum of quasiholes (spectrum at lower filling) by a gap separated into two manifolds: a high energy uninteresting 



2 



set of spurious states, and a low-energy set of quasiholes whose counting empirically matches that of the quasihole states 
in an Abelian i/ = 1/3 filling FQH state. Similar results were obtained for the particle entanglement spectrum a 
procedure used to obtain the topological part of the excitation spectrum directly from the groundstate wavefunction 
[T^ . In the FCI, we have empirically found that the counting of states in the lower energy manifold, be they 
groundstates or quasiholes, exhibits a peculiar, yet un-determined structure. On a finite size x Ny lattice, a 
translationally invariant Hamiltonian for electrons has a spectrum classified by total many-body lattice momenta 
T^fj^i^ix ( mod 27r) and X^iii ^i;/ ( mod 2it), with each particle momentum ki^^y taking values 27^^ /Nx, '2,1x1/ Ny, 
j = 1 . . . Nx, I — 1 . . . Ny. We have found that the counting of states in the low energy (or entanglement energy) 
manifold per momentum sector exhibits multiple degeneracies, although the states themselves are not degenerate - 
there are no other exact degeneracies besides the point-group enforced ones, such as inversion. In many cases, we 
have found that an empirical Pauli principle borrowed from the FQH physics can sometimes (but not always) explain 
the degeneracies observed. While the degeneracy of some model wavefunctions has been worked out for the FQHE 
on a lattice [13], these results cannot be applied to the FCI. 

Our prior analysis strongly suggests the existence of an emergent symmetry for FCFs at rational filling factors. In 
the present manuscript we show that such a symmetry is a many-body translational symmetry of the type exactly 
present in the FQH effect. We show analytically that, if the FCI state is in the same phase as the FQH state 
(including in the presence of quasiholes) , the spectrum separates into quasi-degenerate center of mass sectors which 
exhibit identical counting. Within each center of mass sector, the Hamiltonian eigenstates can be further classified as 
eigenstates of an emergent (in the thermodynamic limit) set of many-body relative translation operators. Using the 
recently obtained Girvin-MacDonald-Platzman (GMP) 14] algebra for Chern insulators |15j for the one-body density 
operators in the valence band, we build a set of many-body relative translational operators which diagonalize the 
Hamiltonian eigenstates in the lower manifold in the thermodynamic limit. We then establish an analytic mapping 
between the counting of zero-mode states in the FQH 2-D BZ on the torus (quasiholes) and that of the FCI levels in 
the low-energy manifold. We show that the counting of quasi-degenerate levels below the gap for the FCI should arise 
from a folding of the zero-mode states in the Fractional Quantum Hall system at identical filling factor. In the process, 
we show how to count and separate the quasihole excitations of the Read-Rezayi (RR)[TB] series in the Fractional 
Quantum Hall effect into two-dimensional BZ momentum sectors, and then how to map these into the momentum 
sectors of the FCI. We provide numerical verification for our analytic results by checking the emergent symmetry 
at work for Laughlin 17. ■ Moore- Read (MR) and RRpLSj states on the checkerboard model of a Chern insulator: 
we predict, from the FQH quasihole counting plus our FQH to FCI analytic mapping, the number of non- Abelian 
quasiholes per momentum sector in the FCI, and then check that it matches the numerical data obtained by exact 
diagonalization for large system sizes. Our work then also shows, as a proof of principle, that Non- Abelian FCIs exist. 
The emergence of the low-manifold of states with folded FQH counting can be regarded as both a consequence as well 
as an imprint of the existence of the topological phase. Its existence is rather mysterious, as the Girvin-MacDonald- 
Platznian one-body algebra from which the many-body generators are built is not valid at large momentum [15j and 
is not in one-to-one mapping with the algebra of the continuum FQH problem. 

The paper is organized as follows. In section 1 we analyze the interacting Landau level problem on the torus in the 
continuum and obtain the set of many-body relative and center of mass symmetries that can classify the interacting 
spectrum by two-dimensional BZ oi N x N momenta where N — GCD{N^, Ne). GCD stands for greatest common 
denominator, iV^ is the number of fluxes that pierce the FQH torus due to the presence of the magnetic field, while N^ 
is the number of electrons. This section has no new material, most of the end results having been obtained by Haldane 
in a classic PRL [Tni, but it does derive the physical equations in much greater level of detail. We then re-formulate 
the translational symmetry in terms of density operators, allowing a closer analogy with the FCI. We also explicitly 
show how to implement the many-body symmetries to build the FQH Hilbert space. Section 2 describes the counting 
of the zero modes of several pseudopotential Hamiltonians per two dimensional BZ momentum. We focus on the 
Read-Rezayi series and show how the generalized Pauli principle that allows the counting of the quasihole states in 
the one-dimensional orbital space on the torus can be tuned to give the counting of zero modes in the two-dimensional 
BZ. As such, resolving the zero-modes (quasiholes) of model Haldane pseudopotential Hamiltonians in the 2D BZ 
becomes a combinatorics problem of counting partitions. In Section 3 we first re-derive the commutator algebra of 
the projected densities in the lowest band of the Chern insulator on a iV^; x Ny lattice first derived in 15j which 
reduces, in the long wavelength limit to the GMP algebra for the FQH effect for N^p — N^ ■ Ny. We show that the 
non-commutativity of the projected densities which gives rise to the GMP algebra is a requirement in the Chern 
insulator, while in the trivial insulator the projected densities are adiabatically continuable to commuting variables 
in the atomic limit. We then show that the existence of a GMP algebra superimposed on a lattice of N^ x Ny sites 
implies the existence of many-body relative translational operators which classify many-body states by a 2-dimensional 
momentum in a reduced BZ of GCD{Nx, Ne) x GCD{Ny, N^) momenta. If the FCI phase is identical to that of the 
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FQH, there is an emergent, center of mass degeneracy of Nx/GCD{Nx, Ne) in the x direction and Ny/GCD{Ny, Ng) 
in the y direction. We show that the counting of states in lower energy manifold of the FCI is a folding of the 
counting of zero-modes in the N x N BZ of the FQH times the difference between the center of mass degeneracies. 
We establish the analytic mapping between the FQH and FCI counting, and give several examples of this counting for 
several Abelian and non-Abelian RR states based on the generalize Pauli principle. We refer the reader versed in FQH 
translational symmetries to read this section and look at Figures [2] [sj [i] for a simple understanding of how the FQH- 
FCI map works. In section 4 we engage in extensive numerical calculations of both energy and entanglement spectra 
that show the existence of Laughlin, MR and Z3 RR states on the checkerboard model |2Qj of the Chern insulator. We 
then show that the counting of quasihole states per momentum sector obtained in numerics matched the one derived 
by our analytic map for all our large set of data and parameters (lattice aspect ratios, electron numbers) tried. Up 
to the largest sizes available on today's computers, the numerical data supports our analytic result. 



MANY-BODY SYMMETRIES OF THE INTERACTING ELECTRONS IN A MAGNETIC FIELD 

In this section, we aim to analyze the symmetries of the interacting Hamiltonian of electrons on a two-dimensional 
torus in the presence of a magnetic field and electron-electron interactions: 

with Ilj = —ihVj — eA{rj) = —ihVj + \e\A(rj) the canonical momentum in the presence of a magnetic field. TVg is the 
number of electrons in the system. We choose not to ^auge-fix, and have V y. A = B. The positions of the particles, 
{rj}, reside on a two-dimensional torus of generators Li, L2. The Hamiltonian is periodic under translations by these 
vectors, V{fi — rj) = {fi — Vj -\- Li^2) and can be written as a sum over the allowed reciprocal vectors q: 

E ^(^^' - ^1^) = 2:4 E ^(-^ E ^'''•'-'''^ (2) 

i^j q i<j 

where A ~ \Li x Z2I is the area of the unit cell, and q = mqi + nq2, m,n and qi = x ^, <72 = ^ ^i- 

We now analyze the single and many-body translational symmetries of the above Hamiltonian. 



Guiding center coordinates and translational symmetries of the 1-body problem 

We first consider A^g = 1 problem and try to find the symmetries of the 1-body problem. A perfectly good 
translational operator could be e*° '^ - which would translate the single-body wavefunction by a. However, this 
operator does not commute with the Hamiltonian because 

[n^,Uf,] = -ihF^^ (3) 

(where Fq/j = daA^ — dpA^ is the (magnetic) field strength applied on the sample). Hence the Hamiltonian wave- 
functions cannot acquire quantum numbers under this operator and one finds an operator which commutes with 
Hamiltonian: 

ATa = n„ - ^(z X r)c< (4) 
where I = yJU/eB is the magnetic length. The commutation relations read: 

[ATa, Kp] = 2i — tp^a - ihe{daAp - dpA^,) = i^ezap (5) 

where we have used Eg = ^teapFap as the uniform magnetic field applied on the sample. By direct calculation, we 
have for the commutators: 



(6) 
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K is called guiding center momentum. Hence the operator that implements the magnetic translation is 

r(a)=exp(^a.i?) (7) 

If translations by different vectors would mutually commute, we could form momentum eigenstates of the system. 
However, they do not. Using the operator relation e^^^ ~ e^^e^e^sl^'^l^ valid when [A, [A,B\\ = [S, [A, i?]] = 
(which is the case here as the commutator of two guiding center momenta is a constant for uniform B), we find: 

T{a + h) = r(6)T(o)e-5^"°^^[^°'-^'*l = T{b)T{a)e- ^'-^^''^'^ (8) 

which leads to 

[T(a),r(6)] - -2sin(i^i. (a X 6))r(a + 6) (9) 

This is called the magnetic translation algebra. This algebra leads to the quantization of flux passing through the 
lattice: going around the unit cell must give us the identity (same condition as commutators of Li, L2 translations 
must commute). We now particularize (without loss of generality, since one can always deform the BZ to a rectangle) 
to Li = LxX, L2 = Lyi): 

T{L,x)T{Lyy)T{-L,x)T{-Lyy) = 1 = e^^^^" (10) 

Hence L^Ly = 2ttPN^ or, for a non-rectangular lattice, -pz ■ {Li x L2) = 2ttN^. For this quantization condition, 

we have that [T(Li), T(L2)] = 0. N^j, is the number of magnetic flux quanta passing through a unit cell. As a 
spoiler for the Chern insulator section, notice that this is the same as the number of sites Ng on a square lattice 
with lattice constant oq = VSttP and number of sites in the x direction L^/aQ and in the y direction Ly/ao- Hence 
the magnetic field to Chern insulator lattice analogy is: L^/ V^irP — > N^, Ly/V^nP — > Ny,N^ Ng. Observe that 

— > Ns are all integers. 



Many-Body Translational Symmetries 

We now focus on the interacting problem. The many-body problem is characterized by translation operators for 
each i'th particle, Ti{a). Translational operators of different particles commute: 

[T,(a),T,-(6)] = -25,jT,{a + h) sin(^5 • (a x b)) (11) 
Physical quantities are invariant under magnetic translations 0/ any particle i by a linear combination of multiples of 

Lx; Ly'. 

Lmn = {mLi + nL2, m, n G Z} (12) 

Following Haldane I!) . primitive translations Lmn are those for which XL„im < A < 1 are not lattice vectors - hence 
{m,n) have no common divisors. Choose a potential V{r + La) = V{r) such that the whole Hamiltonian commutes 
with the translation operators Tj(Li), Ti{L2), for each particle i — \ . . . N^. This means that the wavefunctions of the 
Hamiltonian are eigenstates of the translation operator up to a phase: 

%{L,)ij^{n) = ^^{n + Lj) = e'^^^in) (13) 

where i is the particle index, j = 1,2 and 0* are the eigenvalues. In the wavefunction i^ai^i), we have suppressed 
the position of all other particles but the i'th, as they do not get translated by T^. Let Ti{Lmn) be the translation 
operator that translates the i'th particle by the primitive translation Lmn — niLi + nL2: 

T,{Lmn) IV-a) = 6^'"" IV'a) (14) 

The eigenvalue ^^n, which we will encounter again, can be expressed in terms of the eigenvalues of the primitive 
translations by Li , Z2 , 

T,{Lmn) = T,{mLi + nL2) = r,(mLi)T,(nL2)e*"'^^^^^^^ = (r,(Li))"(r,(i2))"e"^*'"" (15) 
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Hence the eigenvalue of translation by Lmn is Ol^„ = TrmnN^ + m9\ + nO^. Since eigenstates are symmetric or 
antisymmetric under exchange of identical particles i, j, we have that e'^™" = e'*"»" — e*^™": the eigenvalues of the 
translation operators for particle i do not depend on the particle chosen. The Hilbert space of the problem is separated 
in different sectors labeled by 0i,92- At this point, the analysis of the symmetries of the problem under translational 
invariance by each particle has run its course. Haldane |19j showed how to now introduce the many-body formalism 
for the translation operators. We expand his description in detail below. 

Consider electrons on a torus pierced by fluxes, such that the filling factor v — N^/N^ is a rational number 
p/q, with (p, g) = 1 relatively prime. Let iVe = pN , = qN where GCD{Nf., N^) = N. One thing we can hope from 
a many-body formulation of the translation operators is to separate the many-body problem into its center of mass 
part and a relative coordinate part. The center of mass operator moves each and every particle by the same amount: 



TcM{a) = Y[T,{a) (16) 



The center of mass translation by an arbitrary quantity commutes with the Hamiltonian but does not commute with 
the single-particle translation operators Tj {Lj ) and hence does not keep us in a Hilbert space specified by 6*1 , 02 . To 
find the center of mass operators act within the same Hilbert space, we look for all a for which the center of mass 
operator TcM(a) commutes with every Ti[Lj). Since operators of different particles commute, this is tantamount to 
imposing the constraint: 

[T,(a),r,(L,)] = -2r,(a + 4-)sin(^i • (a x L,)) = (17) 

Hence the condition for vanishing commutator is £ • (a x Lj) — 2PiTr. Since Li, L2 span the 2-D plane, it is clear that 

we 



17 



any vector in the plane can take the form a = aLi + PL2 where a, /3 are real numbers. Taking j = 2 in Eq 
have: az ■ (Li x L2) = 2TTpr = a2TrPN^. Hence a — r/N^, and similarly - by taking j = 1, for /3. Hence the most 
general center of mass motion which leaves the system invariant is TcMiLmn / N^) . All other forms of center of mass 
translation change the Hilbert space of the system. 

Our purpose now is to separate the motion of the system into a center of mass motion and a relative motion of 
the particles. Towards this, Haldane [H] defined the relative translation operator Ti acting on particle i so that the 
motion of that particle is compensated by the motion of all other particles in the opposite direction: 



f,{a) = l[T,{a/Ne)Tj{-a/Ne)^T,{{N,~l)a/Ne) Y[ T,{~a/N,) (18) 

J = l j = l, 

Since Ti(a)Ti{—a) — 1, the relative translation operator of the system has the property that Hil^i Ti{a) — 1. When 
Ti{a) applied to a function of fi — rj {i ^ j), it translates the function by a; when applied to a function of Vj — fk 
{k, j, ^ j), it does nothing. Due to the periodicity of the interaction term, only Ti{Lmn) commute with the Hamiltonian. 
We have introduced the center of mass operator and the relative translation operator, so it seems natural that the 
total translation operator factorizes in a product of the two: 

We ^ ^ ^ ^ 

The relative translation operator commutes with any center of mass translation for any a,b\\Ti{a),Tcyi{b)\ = 0, a 
clear indication that both of them arc diagonalizable at the same time. To review the bidding, we have introduced 
center of mass translations T'cm(<J) which commute with the Hamiltonian, found the ones (7cM(^mn/A^0)) which 
commute with the single particle translation operators Ti{Lj) - even though, as we will see, the TcM{Lmn/N^) do 
not necessarily commute between themselves and found the relative momentum operators which commute with the 
Hamiltonian and with the center of mass translations. 

We would like to simultaneously diagonalize the relative translation operators and the Hamiltonian, and are hence 
after the maximal set of Ti{a) which commute with each other [Ti{d),Tj{b)] = and which act in the same Hilbert 
space, i.e. which commute with Tj{L„in). What are the d,b that satisfy these equations? To solve these constraints. 
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it is simple to expand a = aLi + which is always possible, with a, (3 real numbers. From the commutator or 
relative translations with the single particle translations of primitive lattice vectors, we have, for j = i: 

= [f,(a),T,(i„,„)] = ulu. ™^.r„,(-|:)[r,((^^),r,(i„,„)] = 

= -2»n^Li,„^.r„(-£)Ti((^^ +L„0sin(2^^f • (a X (20) 

Hence we find: ^^^j^^z ■ (a x Lmn) = 'iiTpr where r is an integer. To find the constraints on a, (3, we can pick 

Lmn = Li,L2 to obtain: ^^^^^ N^a = r , ^^'^'^^ N^j3 = r. As we are at filling Ng/N^ = p/q, we have: {Ne — l)qa = pr, 
{Ng — l)qp = pr where r is any integer. Also from the commutator or relative translations with the single particle 
translations of primitive lattice vectors, for j ^ i, we have: 

0=[fj{a),Ti{Lmn)] = 
~ iVe ^ ) nm=l, m^i,j '^m{ — Jf^)[Ti{—^),Ti{Lmn)] = 

= -2iT,((^^)n;;;e^^_^^. .T„(_^)T,(-^ +i„Jsin( (21) 

The conditions the above equation gives are: qa = pr, ql3 = pr where r is any integer, consistent with, but more 
restrictive than the first set of conditions. We now impose the condition that the = [Ti{a),Tj{b)]. For i = j we 
have: 

f,(a)f,(6) =T,(^a)T,(^6)a'^'ir,(-^a)r,(-^&) = 

= e-^^-("^>^^^)(^+("-^)^V,(^(a + 6l)n;^,T,(-^(a + 6)) = 

= e-ii--(«x^)(^)T,(^(a + 6))n;^^TK-^(5+6)) (22) 

For f,{a)fi(b) to equal fi{h)f,{d), we need that z ■ {a x b){ ^'\/^ ) = 2Trfr. If a = aLi,b = PL2 we have: 
^^^^^ a(3N^ = r and hence {Ng — l)qaP = pr for r some integer. For i ^ j we have: 

mm=T,{^a- ^h)T,{^h- A_5)n,^^^^.T,(-^(« + 5))e-^'-^'^^^»'-^^^=-^^^+^^^-^'^l (23) 

and hence (by taking Tj{b)Ti{a) and requiring the vanishing of the commutator), we must set z ■ {a x b)[—2{Ng — 
1)^ + {Ne-2)^]= 2TTpr If a = aLi , b = /3L2 we have qajS = pr. 

We now bring the four sets of constraints on a, (3 together. The conditions for the set a for which the relative 
translation operators commute between themselves and commute with the Ti{Lmn) arc: qa = pri, qf3 = pr2, qaf3 = 
Wii 7*11 ''2, ''s G Z If g is a prime number, then we can show that a,/3 = p: we substitute a,/3 from the first two 
equations in the third to obtain pr\r2 = qr^. Since p, q are relatively prime then rz = pr and rir2 = qr where r e Z. 
But since q is prime, then rir2 = qr implies ri = qr' where r' G Z (or vicc-vcrsal for r2). Wc then plug this in 
qa = pri to obtain a = pr' and so a is an integer proportional to p. The smallest value (which gives the largest 
set of a = aLi) of this is a = p. Then qp = r^ = pr and hence (3 = p. The set of vectors a for which relative 
translation operators commute with each other and with the single particle momenta is a = pLmn if ? is a prime 
number. In passing, we note that if q is not a prime, other possibilities arise for the set of maximally commuting 
relative translations. The simplest case is to assume q = qf, i.e. 5 is a perfect square, but same situation occurs 
whenever q is not prime. Then if we choose a = /3 = p/qi wc have ri = r2 = qi,r3 = p - all integers and hence the 
relative translation operators commute with themselves and with the one-body translation operators by Lmn- This 
then gives a = ^Lmn- This, however, corresponds to a difi'erent choice of resolving the groundstate center of mass 

degeneracy, as will be now shown. The center of mass operators TcM{Lmn/N^) commute with the the Ti{Lmn) but 
not between themselves (so they cannot all be simultaneously diagonalized): 

TcM{LmnlN,)TcM{Lm'n'/Ns) = I^. ((£,„„ + L^'n' )/^.)e" ^ ^ ""^'^l = 

= e— ^(""'-""') mLmn + Lm'n')/N,) = e--f n.=i mLmn + Lm'n')/Ns) (24) 
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As a hint of the q-io\d degeneracy of the spectrum, for q a prime number, we see that TQ^^{Lmn/Na) commutes with 
^CM (-^mn)- Moreover, if q is not a prime number, and it is, for example, a perfect square q = q^, we see that a set of 
mutually commuting operators is T^^j 

We now try to diagonalize the maximum commuting set of relative momentum operators Ti{pLmn) which commute 
with themselves and also with the T.i{L,„n) operators. We would like to find its eigenvalues, labeled by a 2-momentum 
k. Following Haldane, we make the assumption that the many-body state (which is an eigenstate of Ti{pLjnn)) 
experiences, when acted on by J2i=i exp(iQ • fj), an increase in its momentum by Q, as long as Q is a reciprocal 



lattice momentum exp{iQ ■ Lmn) = 1- In other words, let the eigenstates of Ti{pLmn) be ip{k)), with eigenvalue A^: 



(25) 



By applying fi{pL^n) on the J^i e**^'^' V'(^) ) we obtain: 



i^ik + Q))= fi{pLmn) Ei e"^-''' 



= T,(^pL_)n,^.r,(-%)(e*Q-- \m) 



-P^S^Q-r-.)f,{pLran) V(fc) ) = 



(26) 



We have then Aj^^ = A^e 'iv„'3'^™" and hence A^: = De 'jve'^'^™". ]j can only be a phase independent on k. The 
constant D can be found by requiring that the k = state remains invariant under all translations. 



D 



ij{k = 0)) = fiipLi) i,{k = 0)) =fi{-pL2) ij{k = 0)) =fi{pL2-pLi) ^lj{k = 0) 



(27) 



However, the last term fi{pL2 — pL{) can be re-expressed in terms of the product of the translation operators of pL^ 
and pLi. We have 

f;(p£2)Ti(-pii) =Ti(^pL2)Ti(-^pLi)n;^iT/(-^pL2)TK^p£i) = 



TVe - 1 „2 s 



Ti{pL2 - pLi)e2 



p'i-^-i'N,) - fi{pL2-pL,)e^'^P"^-^^'"'''^ = TAnL. - V.-^^-Dw- 



12" N. 



Ti{pL2 -pLi)e 



(28) 



This yields the equation D ij;{k = 0)^ = £)2gi(JVe-i)P97r ^ g)^ which gives D = (-l)(^e-i)pg. After establishing 

the constant D, we would like to now find the possible values that k can take. We remark that A^e = pN and = qN 
and hence: 



fi{pL„ 



T.(^£„„)nTK-^) 



(29) 



We know that Ti{Lmn) has fixed eigenvalue 0„,n (which for identical particles does not depend on i). This means 
that the A'''th power of the relative Ti^pL^n) operator has a fixed eigenvalue: 



N, 

.(iv.-i)-Qg 



1 = u^e-*'^-^" 



(30) 



Having determined the constant D, we have = (— l)9^e(JVe i) _ ^ Hence e '^'''^mn = x, and we have as solutions 
k ■ Li = 2m, k ■ L2 = 2'jTj, i,j G Z. We now must ask how many of these solutions represent unique eigenvalues 
of Ti{pLmn)- Its eigenvalues are (— 1)^'*^^="^^ exp(— ifc • Lm„/iV), and are different for k ■ Li = 27ri, k ■ L2 = 27rj, 
i,j G [1 . . . A^]. Hence the BZ of the relative translation operators is made up of Af x A^ values of the momentum. An 
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alternative way of presenting the resulting momenta is to say that the eigenvalues of the relative momentum Ti{pLmn) 
are e~*^'^™"/^, with kl = \J^j;is — sq, A(t — Iq)) where s,t = 1 . . . N and A — L^/Ly is the aspect ratio, while the 

so,to are the quantum numbers belonging to zero momentum: exp(27riso/-/V) = exp(27rto/A^) = (_i)p9(^e-i). ^pj-^jg 
exhausts our discussion of the relative translation operators. 

Having fixed 9i,62,k, we now ask if there are any other degeneracies? The answer is yes. Physically, this is because 
9i , 02 define the single-particle Hilbert space while k is the principal symmetry quantum number of the relative 
wavefunction \iprei)- Left untouched so far is the center of mass translational symmetry of the problem, or the center 
of mass wavefunction \4'cm)- In principle, we have diagonalized only the Ti{pLmn) but there are missing center of mass 
operators, TcMiLmn/Ns) which also commute with the Ti{Lmn) operators and keep us in the same Hilbert space; also, 
any center of mass translation operator commutes with the relative translation operator, and hence we are well on our 
way towards finding other commuting operators. However, crucially, two center of mass translations TcMiLmn/Ns) 
and TcM{Lm'n' /Ng) do not commute with each other and hence they cannot be simultaneously diagonalized. We 
hence must find the maximum set of TcMiLmn/Ng) that can be diagonalized, which, per the above is equal to the 
maximally commuting set of TcM{Lmn/Ns)- We can find out how many TcM{Lmn/Ns) are self-commuting by either 
brute-force calculation or by a smart argument. We start with the brute-force calculation of the commutator: 

[TcM(^),TcM(^)]=e - ||T,( ) = e ^ [[%{ )(31) 

i—l 2 — 1 

which vanishes iff {mn' — nm')/q G Z. Hence there are q values possible. If we require IV'ca/) to be an eigenstate of 
Tcm{L'^ /Ns) where L° is some particular primitive translation, the set of all center of mass translations that commute 
with TcuiL'^ /Ng) is given by {TcM{{qL„in + rL^)/Ns)} with r = 0, 1, ... (7 — 1. Indeed, the commutator 

[TcM{{qL^n + rL°)/iV,),TcM((gim'„' + rL'')/N,)] = (32) 

As Haldane mentions 19 , this is only one of the few resolutions of the center of mass degeneracy. The smart and 
quick argument which reveals the q-fold degeneracy is the following: we know that 

Tj(pi„i„) ^TcM{pL,nn/Ne)fi{pL„in) (33) 

and hence the eigenvalue of 

(TcMiLmn/Ns))'' = TcuipLmn / N^) = Ti{pLmn)T^{-'pLmn) (34) 

is fixed once we have diagonahzed 9i,92,k to be {—iy^'^^'~^^ exp{ip9mn) exp{iqk ■ Lmn/Ns). Since only the g'th 
power of the eigenvalue is fixed, we see that there must be q center of mass operators (of different eigenvalues). If 
the eigenvalue of Tcm{L° /Ng) is A, then the eigenvalues of the set of maximally commuting center of mass operators 
{TcMiiqLmn+rL°)/Ng)} with r = 0, 1, ... g- 1 have eigenvalues X^(-l)Pi^iN^'^) exp(ip0,„„) exp{iqk- Lmn/N,). This 
completes the many-body theoretical symmetry analysis of the spectrum of Fractional Quantum Hall states. The 
section so far did not contain new material, although we believe and hope that the detailed and expanded description 
of the calculations present in [19] is useful to the reader. 

Building the Hilbert space 

On the torus, each Landau level has an identical number of states that it can accommodate. The operator for an 
electron at position {x,y) is ^{x,y) — j 4'm,j{x,y)cmj where Cm,j is the annihilation operator of an electron of 
momentum j in the m'th Landau level, and the (jfm.j are the single-particle orbitals: 

fcez ^ 

Hm is the Hermite polynomial. In the above, we have picked the Landau gauge A = Bxy. While the single-particle 
orbitals given above are not normalized, it is crucial to notice that the normalization factor depends only on the 
Landau level index m and not on the orbital momentum j. This is a feature of the torus geometry, on the sphere 
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the single particle orbitals depend on the angular momentum quantum number. For the lowest Landau level, to 
which we particularize, Hq = 1. The translational properties of the single-particle orbitals (for a rectangular lattice 
Li = L,j.x,L2 = Lyij) are trivially obtained: 

4>m..j{x,y + Ly) ^ (t>^^j{x,y), ^^^(a; + L^, y) = e*"^^*^(/)™j(a;, y) (36) 

In order to obtain the action of the relative translation operators on the many-body states, we note the translational 
properties under translations by {N^ — l)pLmn/N(, and by —pLmn/Ne- 

'Pm,3{x,y+ ^ pLy) = (f),n,]{x,y- -^pnLy) = e^'^^<j),nj{x,y) (37) 



brn,j{x + ^ pLx,y) = e'-y'^'^'^' ^^^0„j+g(x, y), (l)^,j{x - ^pL^,y) = e ' "-y"^" (j)jn,j+q{x,y) (38) 



In the Landau gauge used, the guiding center translation operator T{L) is related to the usual trans- 
lation operator t{L) — expL-V by T{L) = exp{-p{Lxy + \LxLy))t{L). On a one-body state 
T{L)\m,j) = J J dxdyT{L)\x,y) {x,y\m,j) = J J dxdy exp{f2{L^{y + Ly) + ^L.j;Ly) \x + L^,y + Ly)){x,y \m, j) 
where {x,y\m,j) = (j)m.j{x,y). By switching variables in the integral, we h ave : T{L)\m,j) 



J J dxdy exp{-p{Lxy + ^L^Ly) \x,y)){x — L^Ty — Ly \m,j). Using the properties of Eq 37 38 , we can prove the 
following: 

f^ [pLyy) I ji , . . . , ) = e^^ J. | , . . . , j^^ ) (39) 

This was possible because there is no in the L = pLyfj translation, and hence the factor exp(j|-(La;?/ + ^L^Ly) — 1. 
The relative translation operator by pL^x, TilpL^x) has the factor exp{^ pL^y) — exp{ij^q{Ne — l)yi) if it 
acts on the particle i in the single-particle decomposition, which cancels the factor exp{~ij^q{Ne — l)2/i) present 
in {xi — ^"j^^ pLx^ yi |m, j) ■ Similarly, TiipL^x) there contains operators Tj{—pLxx/Ne). These give a factor 
exp{—Y2-^pLxyj) = exp(— i j^gj/j) which cancel the factor exp{ij^qyj) arising from {xj — ^pLx,yj \ni,j). We 
hence obtain 

fi (jpL^x) |ji , . . . , JN, ) = I ji + q, . . . , JN, + q) (40) 

We hence found that (in the Landau Gauge), the many-body Hilbert space vectors . . . , JatJ is an eigenstate of 
Ti(pLyy) with eigenvalue e*T^ 5Zi=i dependent on the total momentum 'Yl!i=i Ji( niod N). Due to the denominator 
N , the momentum is defined only mod N , which is an explicit way of seeing that the relative momentum BZ is 
made out of ky momenta. Note however, and this is essential, that this does not imply that all states with identical 
Si^i Ji( ^) belong to the same Hilbert space. In fact, only states with all Y^^=i3i{ ^<t>) belong to the 

same Hilbert space. The construction of relative translational symmetry-sorted Hilbert space proceeds as follows: 
First, write down all possible states ■ ■ ■ , ,jN^) with no constraint other than no double occupancy of orbitals 
in the case of fermions. Now sort the Hilbert space into different sectors given by the constraint that each sector 
contains terms with identical X^il^i Ji( -^0)- two states having different '}2f=i3ii. '^od A^^) can be coupled 

by a momentum-conserving Hamiltonian. Since — qN , there will, in general, be more several sectors, with different 
Siii ji{ mod Nij,), which have the same X^ill^i Jii mod N). It is hence a mistake to first sort the Hilbert space first by 
sectors with identical X^i^i Ji( niod N): this will result in a much larger number of free-many-body states per sector, 
because the states have yet to be sorted by X^i^i Ji( mod N^). Once we have sorted the states by X^i^i Ji( mod Nfj,), 
for the elements in one sector (which has a momentum exp{—ikyLy/N) = e.xp{i2n^^J:^ii/N)exp{mpq{Nc. — 1))) we 
must implement the translational symmetry in the other direction. As the translation operator TiipL^x) takes 
I ji , ■ • ■ , 3n, ) into I ji -I- g, . . . , jAT^ -I- q) (which crucially, has a the identical Ya=i Ui + q) ( mod N^) = {qNe + Y^f^i 
mod N^) = X^iii jii mod N^), and same thing with replaced by N, so ky remains unchanged), we must form the 
orbits of this operator. This operation goes as follows: for every element of the set . . . , jatJ form all the elements 
\ji + qk, . . . tJNc + 1^) with fc S Z. These elements form the orbit of Ti{pLxx). Let the number of elements in an 
orbit be Z. We can form eigenstates of the Ti{pLxx) by taking combinations of the states in an orbit: 

z 

\kx,ky) =Y,(i'^^\n+(lk,...,jN^+qk) (41) 

fc=0 
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where n is an integer. The ky momentum of these states has been already estabhshed by the equation 
exp{—ikyLy/N) — cxp{i2Tr'^^J:-^ji/N)cxp{iTTpq{Ne — 1)). To estabhsh the kx momentum, we note that 

f^ipLxx) \kx,ky) = e'"-^ \kx,ky) = (_l)P«(^-i) exp{-ikxLx/N) \kx,ky) (42) 

and hence the exp{—ikxLx/N) = exp{i27Tn/ N) exp{iTrpq{Ne — 1)), where n £ [0, . . . , — 1] is an integer. These 
are the basis states exphcitly translationally invariant with the relative translation operator. They allow a large 
simplification of numerical computation and explain the degeneracies observed in the spectrum of any translationally 
invariant Hamiltonian. We will come back in more detail to the explicit construction of a translationally invariant 
Hilbert space in the section [] 

The Density Algebra in the Lowest Landau Level 

The results above can be put in an equivalent form when written in terms of the projected density algebra. In a 
similar way to how one defined the guiding center momentum — Ha — -p{z x f)a, a guiding center coordinate Rq. 
(one-body operator) can be defined: 

Ra = ra- -^(nx z)a, R = -z X K (43) 

where r is the position of the particle. The guiding center coordinate is then related to the guiding center position 
and one can obtain the similar commutations: [Ilci?^] = 0, [Ro^Rp] = il'^Capz- We can not build a guiding center 
density operator: 

p{q) = e^«--^ (44) 
It is known that this operator has several interesting properties. First, it satisfies the GMP algebra: 

p{<l)p{cf) = P{q + <z')e-5''(^~><«"')-, (45) 
which is not surprising as it is just a re-statement of the guiding center translation operator: 

T{L)=pC-^) (46) 

The relation between the translation operators and densities means that all the many-body translation operators 
obtained in the previous section can be expressed in terms of the density operators. If we define the two fundamental 
lattice momenta qi ~ l-nLij A^i^i — —2'kLi/A, where A is the area of the BZ A — 2ttPN^, we obtain r(L,„„) = 
{p{qnm.))^ with q„in = fnqi + nq2, and [(/^((fnm))^* , -ff] = 0. From here on, one can just copy the many-body 
formulation of the translation operators to obtain it in terms of the guiding center densities, and define a center of 
mass density p{qnm) = YIaJLi PiiQnm) and a set of self-commuting relative density operators Pi{pN^qnm) = Pi{{Ne — 
i)pN^qnm/^e)Y[f:^i Pji~P^<l>'lnm/Ne) which Can bc simultaneously diagonalized to obtain the relative momentum 
quantum numbers, an identical copy of the situation encountered in the translation operator example. 

It is crucial to understand how the guiding center density relates to the usual density operator po{q) — exp{iq ■ r). 
It is easy to prove that, up to a factor, the guiding center density equals the projection of the usual density operator 
into the lowest landau level: 

Pe'«-"> = e-^p(g) (47) 

In other words, the guiding center density is the projection of the usual density into the lowest band of the Landau 
level insulator. Let us re-capitulate the bidding. The many-body translational symmetries can be expressed in terms 
of a set of many-body density operators that obey a GMP algebra and who are the projection of the regular density 
operators in the lowest band of the insulator. 
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COUNTING OF ZERO MODES FOR THE READ-REZAYI AND OTHER MODEL STATES (AND 
QUASIHOLES) PER MOMENTUM SECTOR ON THE TORUS 

In this section, we present a procedure which allows for the computation of the number of zero modes (groundstates 
and quasiholes) of many model FQH states resolved per each kx,ky momentum sector of the relative translational 
operators TiipLxx), Ti{pLyy). We aim to obtain this counting without diagonalizing the Hamiltonians for which these 
model states are exact zero modes. We aim to give a combinatoric procedure for counting the states that is valid 
and can be performed for any number of particles. All states whose quasihole counting satisfied a generalized Pauli 
principle in orbital space are amenable to our construction. In a future section, we will use this counting and the one 
numerically obtained from the FCI to check the emergent symmetry proposed in the FCI. 

We work with a system of orbitals. If that system was the sphere manifold, the problem of counting the zero 
modes of certain pseudopotential Hamiltonians can be related to the problem of counting partitions satisfying certain 
conditions. The same conditions, subject to another periodicity constraint, are valid on the torus. We now state 
and justify these conditions, which have been obtained in previous studies |21H29j . Fermionic/bosonic many-body 
wavefunctions of Nf, particles on any manifold can be expressed as linear combinations of Fock states in the occupation 
number basis of the single-particle orbitals. Each Fock state |A) can be labeled either by the list of occupied orbitals, 
A, or by the occupation number configuration, n(A). A = [Ai, A2, . . . , A^rJ is an ordered partition of (the z angular 
momentum on the sphere and) the one-dimensional momentum K = J2i=i ^ii mod N^jy) on the torus into Ne 
parts and each orbital with index Xj is occupied in the Fock state. We number the orbitals on the torus from 
to — 1, and hence each A^ belongs to this interval. By definition, A^ > Xj ii i < j, and re-ordering of the basis 
acquires the appropriate number of negative signs depending on whether the state is bosonic or fcrmionic. n{X) is 
the occupation number configuration. It is defined as n{X) = {nj{X),j — 0, . . . N^}, where '^•^(A) is the occupation 
number of the single-particle orbital with (angular) 1 — D momentum j. In the unnormalized basis in the lowest 
Landau level on the torus, 

(zi, . . . , zjvl A) = S'[0o,Ai (zi) ■ ■■■ ■ (f>o.\N (zn)] (48) 

where S is the process of symmetrization/anti-symmetrization over all indices i,j such that A^ ^ Xj. In this section, we 
will repeatedly use special kinds of partition — the so-called {k, r)-admissible partition. For bosons, a (fc, r)-admissible 
partition labels a configuration whose occupation configuration has no more than k particles in r consecutive orbitals. 
For fermions, a (fc, r)-admissible partition labels a configuration whose occupation number has no more than k particles 
in fc -t- r consecutive orbitals and not more than two particles in two consecutive orbitals. For example, [2, 2, 0, 0] is a 
(2, 2)-admissible bosonic partition on the torus of N^j, = 4, while 2, 1, 0, is not. On the torus, due to its periodicity (the 
orbital of ID momentum is identified with the one of momentum N^) [2, 2, 0, 0] is not a (2, 2)-admissible partition 
on the torus of A^^ = 3, as depicted in Fig. [ijFor obvious reasons, we call these admissibility rules generalized Pauli 
principles. Partitions satisfying Pauli generalized principles as the (fc, r) one above, or other more complex ones play a 
prominent role in the physics of the FQH as it was proved that they count the Hilbert space of the quasiholes of many 
model FQH liquids. The number of (fc, 2)-admissible partitions count the number of quasiholes of the Zj. RR state, 
while the generic (fc,r) case counts the number of Jack polynomials states. On the sphere, other concepts such as 
squeezing, etc are important, but they are not on the torus and we will not discuss them further, referring the reader 
to the existing literature [30] On the sphere and on the torus, the bosonic RR sequence at filling v — k/2 here (see 
Sec.[] for other model states), are the unique, highest density zero-mode wavefunctions of (fc-l-l)-body pseudopotential 
Hamiltonians [3T]. Recent workf^I] has shown the RR bosonic wavefunctions -0 to be Jack polynomials J"^ indexed by 
a parameter a = — (fc-|- 1) and the densest-possible (fc, 2)-admissible 'root' configuration [32] n(Ao) = {fcOfcO/cO . . . /cOfc}. 
On the torus, no extension of the Jack polynomials exist, but the correspondence between zero mode and partition 
counting is still valid. One can imagine it does, as the admissible partitions implement a generalized Pauli exclusion 
statistics, which should be independent of the manifold on which the state resides. On the torus, at filling k/2, 
Nff, — 2N^/k for the groundstate and bosonic (A;, 2)-admissible partitions, we obtain the correct degeneracy [k + 1 
for the Zfe RR groundstates on the torus) by just counting the occupation number configurations with no more than 
A:-particles in 2 consecutive orbitals [A: k 0.... k 0], [fc-1 1 fc-1 1.... k-1 1], [fc-2 2 k-2 2.... k-2 2], 
and so on till [0 k k.... k]. Quasihole counting proceeds identically. Counting states through partitions in this 
way is usually referred to as the "thin" -torus limit, although we stress that this counting also applies for the sphere 
(where it can be analytically proved through the Jack-polynomial mapping), as well as for the finite size torus. 

The counting of quasiholes for the Jack polynomial states proceeds as follows. First, in the ID momentum on 
the torus, write down all the (A:, r) partitions (r = 2 gives the RR series) A = [Ai, . . . , Aat^] possible on the torus 
(remembering that the torus is periodic). Their number is the total number of the quasihole states for that flux. 
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TABLE L Admissible (1, 3) configurations (not more tiian 1 particle in 3 consecutive orbitals) for = 2 particles in A^^ = 12 
orbitals. The total number of such configurations (42) equals the number of zero modes of the Haldane pseudopotential 
Hamiltonian whose highest density groundstate is the Laughlin state. We have N = 2 and hence two ky possible momenta, 
0,7r/Lj, 
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a) [2,2,0,0] N^= 4 b) [2,2,1,0] N^= 4 c) [2,2,0,0] N^= 3 

FIG. 1: A schematic view of the partitions on the torus geometry. The shaded regions are the orbitals (4 for figures a and b, 3 
for figure c). Among the three examples, only the situation a corresponds to a (2, 2)-admissible partition on the torus. On the 
sphere geometry both a and c would be (2, 2)-admissible partitions. 

We now aim to resolve the partitions by assuming they represent states and check their properties under the relative 
translation operators. It is trivial to first resolve ky. To obtain all the states (unresolved in kx) at a certain ky we sort 
the partitions based on the values of Y^fjii Xi{ mod N^). Partitions with the same X^il^i ^<t>) belong to the 

same fey, whose value given by exp{—ikyLy/N) — exp(?27r ^i/^) exp{i7Tpq{Ne — 1)). Call such a set of partitions 
{A'^"}. We now must implement the translational symmetry in the x direction with momentum kx, which we do as in 
the previous section. From all the partitions {A'^f} we again form the orbits by taking all partitions [Ai +qs, . . . , Ajv^ + 
qs], forming the orbit of number of elements Z and form the states X]f=o exp(z27rns/iV)[Ai + qs, . . . , Aat^ + Qs] with 
momentum exp(—ikxLx/N) — exp{i2Trn/N) exp{iTTpq{Ne — 1)). Several important things must be taken into account: 
in the elements of one orbit, [Ai + qs, . . . , Ajv^ + qs] , one needs to bring them to the same form of decreasing elements. 
Note that this is not trivial as A; + qs is defined mod N^, and that it could involve sign changes for fermions. Also, 
for some n's the linear combination of the orbit's partitions might not exist, as some states make up single orbits. 
We will see this in the example below. For other orbits, the number of elements might be smaller than N, in which 
case again not all kx momenta can be present. By building the states \kx, ky) one can clearly see which momenta are 
missing. After momentum resolving, we can easily count the number of states in each kx, ky sector. Note that while 
this procedure is constructive (we do not have a theoretical expression for the number of quasiholes per sector, but 
only a way to construct them), it is an analytic map from a series of admissible partitions in II? to a two-dimensional 
relative momentum BZ. 

A complete understanding of how the counting of quasiholes per sector works comes from working out an example. 
We will try to resolve the counting of quasihole states of a = 1/3 Laughlin state of two electrons iVg = 2 on a torus of 
= 12. The quasihole states are the zero-modes of the Laughlin quasiholes. Per our analysis above, the spectrum of 
any translationally invariant Hamiltonian has a (note Ne/N^ = 1/6 is the filling of the system; the degeneracy comes 
from it, but we are looking at the quasiholes of the z/ = 1/3 Laughlin state, and the counting will reflect that) six-fold 
degeneracy. We will aim to show that the counting of zero-modes satisfies this center-of-mass degeneracy. Note that 
the value of Ne/N^ fixes the degeneracy while the specified FQH state whose quasiholes we try to find = 1/3 fixes 
the Pauli principle responsible for the counting. To start, wc write down all the (1, 3)-admissible partitions in N^f, = 12 
orbitals of iVe = 2 particles. There are 42 such configurations and they are given in Table [l]. Since N^/N^ — 1/6, by 
our theory, the spectrum should be 6-fold center of mass degenerate. The non-degenerate states should be resolved 
by ky and kx relative translation momenta, which each take GCD{Ne, N^) — 2 values 0, 2tt /Ly and 0, 2tt/Lx- Per our 
procedure, we first sort the state into ones with values X^i^i niod N^f,). Table [ll] presents the partitions sorted in 
such way (they can obviously also be read from the 5'th column of Table [l] 

Notice that the mod = sector (which is the equivalent of the fc^ = sector, as (— l)P9(JVe-i) _ j-^^^g 3 

elements while the '}2ix ™^od = 1 sector (the ky — 2n/Ly sector) has 4 elements. Notice how the counting of each 
of the ky = 0,TT/Ly sectors is g = 6-fold degenerate, as the theoretical symmetry analysis mentions. We now would 
like to implement the x translational symmetry to find the counting of Laughlin quasiholes per momentum sector. 
The same result will be obtained no matter which of the 6-fold degenerate sectors we pick. 

We first look at fej, = sector with the states mapped by the partitions [8,4], [9,3], [10,2]. To form the orbits 
we must look at all elements [Ai + qs, A2 + qs] — \\\ + 6s, A2 + 6s]. The orbit of [8,4] contains [14, 10] = [2, 10] 
(and obviously the orbit of [10,2] contains [16,8] = [4,8]). Hence the two states [8,4], [10,2] can be used to build 
momentum eigenstates: 

\n,ky = 0) = [8,4] +e*^ [14, 10] = [8, 4] -f e*^ [2, 10] = [8, 4] - e""[10, 2] (49) 
where we have used N ~ 2 and [2,10] — —[10,2] because the particles are fermions. The action of the relative 
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^ Ai mod N4, 


[Ai,A2] 


E Ai mod N 





[8,4]; [9,3]; [10,2] 





1 


[8, 5]; [9, 4]; [10, 3]; [11, 2] 


1 


2 


[9,5]; [10,4]; [11,3] 





3 


[3,0]; [9, 6]; [10, 5]; [11, 4] 


1 


4 


[4,0]; [10,6]; [11,5] 





5 


[5,0]; [4, 1]; [10, 7]; [11, 6] 


1 


6 


[6,0]; [5, 1];[11,7] 





7 


[7,0]; [5, 2]; [11, 8]; [6,1] 


1 


8 


[8,0]; [6,2]; [7, 1] 





9 


[6, 3]; [7, 2]; [8,1]; [9,0] 


1 


10 


[7, 3]; [8, 2]; [9,1] 





11 


[7, 4]; [8, 3); [9, 2]; [10, 1] 


1 



TABLE II: Partitions [Ai, A2] sorted by identical Ai mod N^. 



translation operator gives: 

fi{pL^x) |0, ky) = [14, 10] - [16, 8] = [2, 10] - [4, 8] = |0, ky) (50) 

and hence this state has momentum = 0. Similarly, the state ll,^;^,) has momentum kx = 21: /Lx- Out of the 
original 3 states at momentum ky = 0, wc now have to analyze [9, 3]. As it is a single state, it forms an orbit by itself. 
This is indeed true, as [9 + 6, 3 + 6] — [15, 9] = [3, 9] = —[9, 3]. Due to the fermionic sign, this is state at momentum 
kx — 2tt/Lx- We have now determined that wo have 1 zero mode state at momentum kx,ky = 0, and 2 at momentum 
kx,ky = 2'!r/Lx, 27T/Ly. Note that of course, our choice of which states to pick for the ky = momentum is irrelevant, 
giving the same counting. If instead we had picked the set of ky = states [9, 5], [10, 4], [11, 3], then [9, 5] — [11, 3] would 
be the state at kx = 0, while [9,5] + [11,3] and [10,4] would be the state at momentum kx = 2t^/Lx. This completes 
the symmetry analysis of the zero modes of 2 particles in A^^ = 12 of the Haldane pseudopotential Hamiltonian for 
the Laughlin state at ky = 0. 

We now look at ky = 27r/i sector with the states denoted by the partitions [8, 5], [9, 4], [10, 3], [11, 2]. The orbit of 
[8,5] is [2,11] = -[11,2], while the orbit of [9,4] is [3,10] = -[10,3]. Hence the state [8,5] - [11,2] has momentum 
kx = Q while [8,5] + [11,2] has kx = 2-k/L^. Identically, [9,4] - [10,3] and[9,4] + [10,3] have momenta k^ = Q^i^/L^. 
Hence there are 2 states for kx = Q and 2 states for kx = 2'k/Lx- The same counting would be obtained if we 
picked any of the other states at ky = 271 /Ly. Wc have now given a combinatoric prescription for counting the 
states per momentum sector of the relative translation operators. For 2-particles in 12 orbitals the quasiholes of the 
Laughlin state have the following quantum numbers: 2 states in each of the kx,ky equal to 0,2TT/Ly, 2tt/Lx,0, and 
2'k/Lx. 2ii/Ly momenta, and 1 state in the 0, momentum. There is 6-fold center of mass degeneracy which multiplies 
these 7 states to obtain the total of 42 quasihole states. 

As a last simple example, let us look at resolving the quantum numbers of the bosonic groundstates of the Haldane 
pseudopotential for the MR interaction for Nf,= A particles in = 4 fluxes. As Ne/N^ = 1 there is no center of mass 
degeneracy and the BZ of the relative translation operators is a 4x4 momenta. We first write down all the {k, r) = (2, 2) 
partitions (whose counting corresponds to the bosonic MR state) of 4 particles in 4 orbitals. They are, in occupation 
number 2020, 0202, 1111, totaling 3, the correct non-Abelian degeneracy of the MR state, and they correspond to the 
partitions A = [2, 2, 0, 0], [3, 3, 1, 1], [3, 2, 1, 0] respectively. We first sort the states by 'Ya=i ^ii niod 4), which gives 
0, 0, 2. The J2t=i mod 4) = 2 state is [3, 2, 1, 0] has e-'*=f = ei2n-2/4 . ^^le - sign coming from pq{Ne - 1) 

being odd. Hence ky = 0. We now find its kx'- fi{pLxx)[3, 2, 1, 0] = [4, 3, 2, 1] = [0, 3, 2, 1] = [3, 2, 1, 0] by virtue of the 
bosonic nature of the state. Hence 6"**^^^="/^ — ~1 or kx = in/Lx- There is 1 state at kx,ky = in/Lx-O- The two 
other states [2, 2, 0, 0], [3, 3, 1, 1] belong to the ky — An/Ly sector and to the same orbit. We hence form the states: 

\0,ky) = [2, 2, 0,0] + [3, 3, 1,1], \l,ky) = [2,2,0,0]+e*t[3,3,l,l] = [2, 2, 0, 0] + z[3, 3, 1, 1] 

[2, ky) = [2, 2, 0, 0] - [3, 3, 1, 1], [3, ky) = [2, 2, 0, 0] + e*^ [3, 3, 1, 1] = [2, 2, 0, 0] - z[3, 3, 1, 1] (51) 

It is clear that two states [2, 2, 0, 0], [3, 3, 1, 1] cannot be made into the 4 independent momentum states above. Indeed, 
this is because [1, ky) , [3, ky) are not momentum eigenstates; for example Ti{pLxx) [1, ky) = i [3, ky). The only two 
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momentum eigenstates are |0, fcj,) at momentum kx = A-k/Lx and |2,fcj,), at momentum kx = 27r(4 mod N)/Lx = 0. 
We now found that our combinatoric counting principle predicts that the zero modes of the Haldane pseudopotential 
for which the MR wavefunction is the groundstate are resolved by the relative translation operators into momenta 
kx,ky equal to {Air/ Lx,0), (0,47r/Ly) and Air /Ix,^!^ / Ly). An exact diagonalization of the Hamiltonian confirms this. 

It is immediately obvious that the current example can be easily generalized to all the {k, r) states. Our procedure 
allows for a combinatoric approach to resolving the number of quasiholes per momentum sector just by counting 
partitions. It goes beyond the usual counting on the ID-torus momentum which can only give the total number of 
quasiholes per 1 D momentum. More examples and data for Abelian and non-Abelian states will be presented in the 
section which uncovers the mapping of FQH to FCI. 

EMERGENT MANY-BODY SYMMETRIES IN THE FRACTIONAL CHERN INSULATOR 

In this section we show the appearance of many-body translational symmetry operators in the FCI. We first re- 
derive the result, obtained in [TS] and rederived in [33], that the long- wavelength algebra of the projected density 
operators in a Chern insulator has a GMP form when the Berry curvature variation is not large. We show that the 
non-commutativity of the projected density matrices is required by the nonzero Chern number in the nontrivial Chern 
insulator. We then show that the problem of a Hamiltonian with such a projected density algebra can be mapped into 
a FQH problem of a translationally invariant Hamiltonian superimposed on a background lattice of Nx x Ny{= N^) 
sites. We then show that this problem admits relative translational operators and momenta which are a folding of 
the ones obtained in the FQH case in the continuum. 

One-body projected density algebra of an insulator 

In this subsection, we re-derive the result obtained in [TS]. We work with translationally invariant one-body 
Hamiltonians: 

^ ^ 4ahi-3;a,l3Cjp (52) 

where a,/? contain orbital and spin indices. We call a an orbital index. We now quickly introduce some conventions. 
Our Fourier transform sign convention is Cka = 'Ylij 6~*'^''cja where by kj we mean k ■ j and Ng is the number of 
sites on a lattice with aspect ratio Nx x Ny. The Bloch Hamiltonian matrix is: h^p (k) = J^r ^ **''^^r;Q.;9 The Bloch 

Hamiltonian can be separated into normal mode operators 7^ at momentum k of the band n H — J^k n ^■ni^)lk^ Ik 
where the normal modes can be written as a matrix rotation of the original electron operators 7^ = u^*pCk,f3- The 
elements of the matrix are u^"^ which form the eigenstates of the Bloch Hamiltonian ^ap{k)u^ p — En{k)u^ ^. 
We then define the projector into the band n at momentum k by 

P„,=7,"t|0)(0|7," (53) 

Fractional topological insulators are usually constructed and observed in models with fractionally filled bands whose 
bandwidth is very small, such that interactions and not the kinetic energy dominate the physics. The ideal example 
of such an insulator is the flat band one-body deformation of any Hamiltonian insulator h(k): 

Hfb = - (54) 

From now on, we will consider the physics of the (fractionally) occupied bands and look only at projectors into the 
occupied bands. We define the density operator per orbital a at site j, = CjaCja whose Fourrier transform reads: 
~ — Sfe k^a,k+q- Noticc that we do not shift the sum to put k + q/2, k~q/2 because, while 

perfectly suitable in the continuum, that would change the boundary conditions on a lattice at which is quantized 
in units of 2n /Ni (i = x,y) 

As we are interested in the low-energy physics in the fractionally filled bands, we define the projected density 
operator Pkl^'k' = Pn.MPqPn^M which has the following form: Pkl'^q.k,' = 7^w"!fe>a'fci+/fe2,fei+97fei'^ |0) (017^'+, 
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The full projected density involves a summation over the orbital a, and over the projection operators P = Pn,k' 

P^-TW ^ <k<:,-,,lT'\0){0\lT+, (55) 
Taking of two densities at momenta q and w, we obtain: 

fc,ni ,n2 ,Tn2 ,a,/3 

We particularize to the low- wavelength continuum limit q << 1 (in lattice constant units). On a discrete lattice, low 
q represents a small integer combination of 2'n/Ni. In this limit, a Taylor expansion gives: 

<,k+q = K.k + QtdiK^k + \<li<l]didjult, (57) 



where d denote discrete momentum derivatives. We expand the difference in Eq 56 in the wavelength limit, and use 
the identities Y.c.<%<m = '^"i."- E ja.<^;)u- ^ + = and 

A^Amul%M<M) + m<%)u7M + <M^d,d,u2,^) = (58) 

(valid for Aij a symmetric tensor), to find, after tedious algebra (presented in Appendix 1), that in the long wavelength 
limit: 

Eni* n2 Ti2* ^2 ni* 712 TI2* ^2 ^ /" \ TjTi-im2 

'^a,k'^a,k+q^l3,k+q^l3,k+w+q ^ P ,k'^ P ,k+w'^a,k+w'^a,k+w+q ~ l^K^i'^j ^ ^iljJ^ij ) 

"2 

where is the non-Abehan field strength: i^"^'"" = diA"^'""^ -9^ A"^'"" -i[Ai, Ajf^^"^^ with the Berry potential 

^nin2 ^ -iY^pufdjU^J^^ the non-Abehan field strength. The reference [15] hence finds: 

[Pq.Pn,] = -l{q^w,-w^q,)^^ E ^"""^7?^ |0) (0|7r+V. (60) 

A:,ni ,m2 

The difference with the result of [15| is only apparent because we are working in the low wavelength limit q, w small 
and we kept terms to second order in q,w^ to the order in perturbation theory that we are working at, we have that 
Tlib '^b*i^+)'^bi^-) ~ Tlib "b*(^+)"b (''') — ^ which proves the equivalence between the formula in 15 and the above. 
The smallest qx,Wy are qx = 2n/Nx and qy — 2TT/wy. 

Before we investigate the full implications of this result, we point out three important results. First, the commutator 
of two densities has to be non-zero in a nontrivial Chern insulator. That is so because the Chern number of the 
insulator can be expressed as a trace over the density commutator, as shown in Appendix 2. Second, the commutator 
of densities at the smallest available lattice momentum 27: / {Nxa),2-K / {NyO) with a the lattice constant is identical 
to the commutator of the projected X and Y lattice position operators, as shown in Appendix 3. The trace of this 
commutator is known to give the Chern number. Third, as a bonus, the expression of the Chern number in terms of 
the commutator of the projected X, Y coordinates allows for a direct derivation of the Chern number of the continuum 
Landau level, which was not possible with the usual TKNN formula for the Chern number which involve derivatives 
over the two momenta. We show this in Appendix 4 



Girvin-MacDonald-Platzman Algebra for Insulators with Smooth Berry Curvature 

We now particularize to two-band models (insulators with one band below and above the gap) or to many-band 
insulators where the non-Abelian components of the field strength can be neglected. We observe that the commutator 



in Eq 60 does not form an algebra as the right-hand side is not easily expressible in terms of a density operator. 
However, the commutator algebra has an eerie similarity to a CMP algebra at long wavelength. As pointed out in 
[15j . if the local Berry curvature can be replaced by its average: 



^^^cPkBmik) _ 2-kC 
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where C is the band Chern number. We then observe that the field strength can be taken out of the sum in Eq 60 , 
and, since to zeroth order m q + w, Pq+w — X]fc 7fc^ |0) (0|7fe+ij+«) + + The commutator algebra then becomes: 

r 1 ^\ - 27rC 

[Pq-, Pw\ = X w) ■ ^ (27r/a)2 ^'^+^ -* 

which is identical to the GMP algebra at small q. Several comments are now necessary: first, in a two-band insulator, 
it is impossible to have a constant Berry curvature due to the no- hair theorem [34j. although this seems possible in 
insulators with four or more bands |34j . However, we will from now on assume that this is indeed possible, and assume 
the validity, at low-energies of the GMP algebra for Chern insulators. We then observe that the density algebra for 
the Chern insulator is identical for low momenta to the one for the Landau level QH problem if: 

27r 

(q^w)- z = {qi X • zl^ (63) 

where q,w are the smallest values of the momentum on the lattice: q = 2tt /(Nxa)x,w = 2Tr/(Nya)y while qi = 
2TTL2/A,q2 = —2'kL\IA are the smallest momenta of the QH torus. A = [Li x L2) • z = 2'kPN^ is the area of 
the continuum torus BZ, (27r/a)^ is the area of the BZ on the lattice, and we specialized, without loss of generality 
to a square lattice. Since the LLL has Chern number 1 (see Appendix 4), we also took C = 1. We hence have 
(gi X 52) ■ zP = {2'kYP /A = 27: /N^ = {q x w) ■ zj^j^ = 2T:/{NxNy) which leads us to the equivalence of the 
GMP algebras for a Chern insulator with Chern number 1 on a lattice with Nx,Ny sites in the x, y directions and the 
algebra of the QH effect on a torus pierced by A^^ flux quanta if: 

N^ = N,- Ny (64) 

An important remark should be made here. As pointed out in |15| . it is important to realize that there is no one- 
to-one mapping between the Chern insulator density algebra and the FQH one. First, the FCI density GMP algebra 
is valid only at long- wavelengths, and it will be our assumption that as system sizes reach thermodynamic limit, 
its application transcends the long-wavelength limit and it becomes an approximate symmetry of the system at all 
wavelengths. This is indeed justified posteriori by our extensive numerical checks, which prove the existence of center 
of mass degeneracies in the Chern insulator which would be possible only if the many-body translational operators 
obeyed GMP algebra for all wavelengths. Second, the number of density generators in the FCI and FQH differ {Ns in 
the FCI, N"^ in the FQH) [15], but, as was pointed out in [15], there are many instances when the significant physics 
is governed only by the operators in the FCI. 



Fractional Quantum Hall to Fractional Chern Insulator Mapping 

We now found that the constant Berry curvature Chern insulator density algebra is identical to that of electrons 
in a uniform magnetic field with number of fluxes = ■ Ny. Hence the translational symmetries of the two 
problems should share similarities. To complete the analogy, since the Chern insulator is defined on a lattice, and 
hence we need to introduce a lattice to the FQH problem. We are then led to consider the translational symmetries 
of the Hamiltonian in magnetic field: 



^ = o!;; E n| + ^ ^ ^ e^?(--~^) + ^^o E E E '^(-1 - - ^^vv) (65) 



j q i<j i—1 m—1 n—1 

with Hj = —ihWj — eA[r-j) — —ihS7j + \e\A{rj) the canonical momentum in the presence of a magnetic field. Ne is the 
number of electrons in the system. We choose not to ^auge-fix, and have W x A = B. The positions of the particles, 
{ri}, reside on a two-dimensional torus of generators Li, L2- The Hamiltonian is periodic under translations by these 
vectors, V{ri — fj) ~ {ri — fj + L1.2) and can be written as a sum over the allowed reciprocal vectors q. A = \Li x L2I 
is the area of the sample, and q = mqi + nq2, ni,n and qi — x z, q2 — —^z x Li. We choose to work on a 



square lattice, without any loss of generality. If the Hamiltonian of Eq 65 was made out of only the first two terms, 
the symmetry analysis would have been identical to the one presented in Section 2 and the number of states of any 
model FQH state per momentum sector would have been obtainable by our constructive procedure of implementing 
the relative translational symmetries in both x^y directions on the (A:, r)-generalized Pauli principle states. 



18 



The last term of Eq 65 implements the lattice with sites in the x direction and Ny sites in the y. The lattice 
constants in the x,y directions are Ix = Li/N^, ly = LijNy It makes electrons energetically favorable (for large 
negative Vb) to sit on lattice sites, and it makes the momenta in the x,y direction take the discretized set of values 
{q ■ x^q ■ jj) — (27ri/ {N^lx), 27rj/ (Nyly)) with i = 0, . . . , N^, j — 0, . . . , Ny, which also serve as the momenta of the 



many-body states. The problem defined by Eq 65 is similar to that of the Chern insulator: it has the same density 
algebra in the long wavelength limit, the same lattice and the same translationally invariant Hamiltonian. The two- 
body Hamiltonian can be replaced by an n-body Hamiltonian, as long as it is translationally invariant. A physical 
assumption is then the fact that the low-energy spectrum of the FQH and FCI problems (when the FCI problem has 
a topologically ordered groundstate) will be similar. More to the point, we take it as a proof of the topological nature 
of the FCI state that the low-energy manifold of states has the identical counting, per momentum sector to that of 



the FQH states at the same filling in the presence of the lattice in Eq 65 . We now find what that counting is. 

The presence of the lattice infiuences the symmetries of the problem. The set of continuum translation operators 
Ti(a) are not valid unless a is a multiple of the lattice constants lxX,lyy. Hence translational symmetry operators 



are Ti{mxlxx),Ti(mylyy). Although the one-body terms of the Eq 65 commute with Ti{lxx),Ti{lyy), due to the 
many-body interacting term (which does not commute with those operators), the Hamiltonian still commutes only 
with Ti{Nxlxx),Ti{Nylyy). The Hilbert space is described by the two twist parameters: 

T^iNJxx) IV) = e'"'' m , TdNylyy) \^) = e^^" |V) (66) 

As the number of fluxes per unit cell plaquette is N^/ {NxNy) = 1, the translational operators of unit cell length com- 
mute: [T,{lxx),Tj{lyy)] = -2S^jT,{lxX + lyy) sin{lxly/2P) = since lxly/2P = LxLy/2PNxNy = 2ttP N^/2f NxNy = 
TT. The center of mass translation operator T(a) = Iliii^il^) commutes with the translation operator Ti(lxx) or 
Tiilyy) if a = IxX or ly-y. Notice that the operator T{Li/N^) (or T{L2/N^)), which commuted with Ti{L2), (or Ti{Li) 
respectively) in the continuum limit, does not appear here because Li/N^ = Li/(NxNy) = Ix/Ny is smaller than 
the permissible lattice constant. Hence the center of mass translation operators commute between themselves and 
theexact g-fold degeneracy apparent in the continuum FQH is reduced to a one-fold degeneracy of any sector (in other 
words, we will not have exact degeneracies anymore). However, if by adding the lattice we still remain in the FQH 
state, then remnants of the original degeneracies and spectrum in the continuum should appear even after adding 
the lattice. In particular, we know the counting of the degenerate zero-mode quasihole states (for FQH series such 
as RR) in the continuum, where they are zero-modes of pseudopotential Hamiltonians. We now try to find what the 
counting on the lattice momentum. 

As we know the counting of the FQH per continuum momentum of the relative translational operators, we now 
try to implement the relative translational symmetries of the problem. The relative translation operators in the 
continuum, Ti{a) = Ti{{Ne — l)d/Ne)Y\^^^Ti{—d/Ne), make sense only if a/Ne is a multiple of Ix or ly, and hence 
the relative translation operators by the smallest amount are: 

f,{NJxx) = T,((iVe - l)lxx)Y[T,{-lxx) (67) 

and similarly for ly. We now ask two questions : do these operators commute with Ti{Li),Ti{L2) and do they 
commute with the Hamiltonian of the problem? The answer to the first question is an unequivocal yes, as Lily/2P ~ 
LxLy/2PNy = nN^/Ny = ttNx). The answer to the last question is no, because Ti{Nf.lxX) commutes with exp(i(f- (r*, — 
ri)) when j,l ^ i as it translates both coordinates by the same amount but does not commute with exp{iq ■ {fi — f))) 

fi{NJxx)e'^^^^-^''> = e*«>"''+(^'^i)'-*)e-*«-(*''-'-*)e'«(^"''-*'')T'i(iVe?^i) = e'^''^^^-^'' e'^^^^~^''>f,{NJxS:) (68) 

since q-x — 211] / Nxlx, j = 0,1,... iVa;, and hence N^lx = 2TTNe/Nx', an identical situation obviously occurs for 
Ti(Nelyy). As Ne, Nx, Ny are integers, we now define several important ratios and common denominators: 

GCD{N„Nx)^Nox, GCD{N„Ny)=Noy, ^ = — , ^ = - (69) 

Nx qx Ny qy 

The qx , qy are integers and should not be mistaken for the components of a g momentum vector. Fortunately, this mis- 
take is hard to do from the context. There are obvious relations between the above and N — GCD{Ne, N^), N^/N^ = 
p/q, but we will come to this later. It is then clear to see which relative translation operators commute with the 
Hamiltonian: 

fi{qxNelxx) = Ti{qx{Ne - l)lxx)Y\_Tj{-qxlxS:), fi{qyNJyy) = Ti{qy{Ne - i)lyy)Y\_^ji^%hy) (^0) 
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These operators commute with both translations and the Hamihonian. They also commute between themselves and 
hence can be diagonalized (remember that in the Chern insulator, the commutations are only approximate, low- 
energy properties, as the GMP algebra is valid only at long wavelength). We now ask how many eigenvalues of these 
operators are independent. We use the same reasoning as in the continuum limit and assume that the many-body 
state experiences, when acted upon by X^^j^ exp(i(5 • rj), an increase in its momentum by Q as long as Q is an allowed 
lattice momentum ex.p{iQxLi) = ex.p{QyL2) = 1. The eigenvalues Afc 

fiiqxNJxX + QyNJyy) iV(fc)) = Xik) iV(fc)) (71) 

can be easily obtained by following the continuum proof closely and satisfy A(fc) = D ex.p{—i{kyqyly — kxQxlx))- 
The coefficient D can be determined by requiring that the A; = state has maximum lattice symmetry, and hence 

D IV^fc = 0)) = f,{qyNJyy) \^{k ^ 0)) = f,{-qxNJxx) \^p{k = 0)) = f,{qyNJyy - q^NJ^x) \^P{k ^ 0)). From here 
we obtain that D = (^—iy^iyNe{Ne-i) _ i (notice the difference from the continuum case, in which D could be both 
1, —1). To obtain the values of different relative momenta, we notice that: 

{f,{qxNJxx)f'"'' = {f,XqyNJyy)f»y = 1 (72) 

which implies exp(— ifc^gx-^Ox^a:) = exp{—ikyqyNoyly) = 1. Expectedly, kx = 2'Kj/Nxlx and ky = 2'Kj/Nyly with 

j = 1,2 Now that we obtained the values of the different momenta, we can ask how many of these eigenvalues 

are distinct. The eigenvalues appear as exp{ikxqxlx) = exp{i2Trjqx/Nx) = exp{i2TTj /Nqx) and similarly for the y 
component, so the number of distinct momenta resolving the relative translation symmetry is: 

kx = P^, i = 0,...,7Vo^-l; ky = ^, I = 0, . . . ,Noy - 1 (73) 
^^x*'x y y 

This defines a Nqx x -^Oy BZ. The momenta of the relative translational operators of the FCI, kx = j = 

0,...,Nox-l; ky = ^, / = 0,...,Afoi,- fare in fact foldings of the fc^ = 1^, j = 0,...,iV; ky = ^, 1 = 
0, . . . ,N momenta of the FQH. To show this, we first note that both Nqx, Noy divide N. Indeed, since GCD{Nx, Ne) = 
Nox, GCD{Ny,N,) = Noy, and iN,,N,) = iNxNy,Ne) = N, we have that GCD{NxNy/Nox, NJNox) = N/N^x = 
GCD{qxNy,pxNox), which proves that N/Nqx is an integer and similarly for N/Noy. With Ne/Ng = p/q, we have 
Ng = pN = PxNqx and hence N = PxN^x/p = PyN^y/p we have that the FCI momentum quantum numbers kx, ky 
of the relative translation operator are a Px/p,Py/p folding of the FQH momentum quantum numbers in the x,y 
directions respectively. We have hence taken care of the relative momentum quantum numbers. 

In the FQH in the continuum there is a center of mass degeneracy of q. This center of mass degeneracy comes from 
the non-commutativity of the continuum center of mass translation operators TcM{Lmn/^s)- On the lattice, we also 
know that since: 

fi{qxNJxx) ■ Tcuiqxlxx) = Ti{qxNJxx) (74) 

the eigenvalue of Tcm{<1xIxx) is fixed once kx and 6x are known (similarly for TcMiQylyV))- Since we have already 
resolved the motion in the Nqx x A'oj, Brillouin zone, we have a remainder of qx x qy momenta to reach the full 
lattice BZ of Nx x Ny. TcM{lxx),TcM{lyy) form a maximally commuting set of center of mass operators which can 
be simultaneously diagonalized (note that in the FQH TcM{Lmn/^s) could not be simultaneously diagonalized) and 
which commute with the Hamiltonian, so no exact degeneracies are generically present. However, if the system is in 
the same universality class as the FQH, the low-energy modes will then exhibit a center of mass qx quasi-degeneracy in 
the X direction and qy quasi-degeneracy in the ky direction. This degeneracy, as well as the relative momentum, make 
up the full BZ of the Nx x Ny lattice. Note that in the lattice example, all translation operators commute with each 
other (courtesy of unit flux per unit cell), so in effect we could just diagonalize those. However, we choose to build 
the relative translational operators because in the continuum case we already know the counting per the momentum 
quantum number attached with these operators in the FQH. To obtain the coimting of the cigcnstates per FCI relative 
momentum, one needs to fold the FQH relative momentum BZ. This folding now needs to be supplemented by the 
miss-match in the center of mass degeneracy. In the FQH, the center of mass degeneracy is q; in the FCI, it is qxQy 
In Appendix 5 we show that q/{qxqy) is an integer number. The counting of states in the low manifold, sorted, per 
FCI relative momentum sector is a folding of the FQH relative momentum BZ times the mismatch in center of mass 
degeneracy q/iqxQy)- 

In conclusion, the diagonalization of a Chern insulator at fractional rational filling should exhibit, in the topologically 
ordered phase, a low-energy manifold of states with counting which encodes the topological character of the state. The 
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FIG. 2: Energy spectrum for TV = 2,Nx = 3, A^'y = 4. The total number of states of the low-energy manifold the dotted line 
is 42. Per momentum sector, the counting is 3 states if ky is divisible by 2 and 4 states otherwise, irrespective of (there 
are degeneracies in the spectrum, for example, there are 5 states at momentum (0,0), the low-energy manifold being 3-fold 
degenerate with zero energy 

low-energy manifold of states is a result of the emergent translational symmetry of the system that the long wavelength 
GMP algebra implies. For a x Ny lattice with electrons, the states in this low manifold will experience an 
emergent center of mass degeneracy, with the spectrum being x qy quasi-degenerate (in the k^, ky direction). If 
the state developing in the system is at filling p/ q, the low-energy manifold will contain either the groundstates or 
the quasihole excitations of this state. Once the center of mass degeneracy is resolved, the spectrum of the FCI per 
relative momentum sector is a folding of the spectrum of the FQH relative momentum BZ. We present clear examples 
of this procedure in the next subsection for the FCI at 1/3- Abelian topological order before presenting numerical 
data on the Abelian and non- Abelian FCI states. 

Two Simple Examples of the FQH to FCI Mapping 

We now look at two simple examples of how the analytical construction in the previous sections can be applied to 
the numerical data in both the energy and the entanglement spectrum. The counting per momentum sector can be 
obtained from our FQH to FCI mapping. Lets first look at the specific case of such counting for N.^. x Ny = 3x4 
lattice Chern insulator model of |20] diagonalized in the flat-band limit [5j with TVg = 2 electrons and nearest neighbor 
(on-site being forbidden) Hubbard interaction. This interaction favors a 1/3 Abelian state, which we know occurs 
in the FQH state for short-range repulsive pseudopotentials. As the state contains only 2-particles, we expect that 
the physics is that of 2-electron = 1/3 Abelian FQH state plus a number of 6 quasiholes. The physics also reveals 
the structure of the pseudopotential two-body interaction in the FCI. The spectrum, shown in Fig[2], confirms this. 
We observe a clear gap between a low-energy manifold of states and a high-energy manifold of spurious states. We 
focus on the low-energy manifold. By looking at the numerical data (Figjs], which includes degeneracies), we see 
that the low-energy manifold exhibits a counting of states of either 3 or 4 states per momentum sector, depending 
on whether ky is divisible or not by 2, as can be seen in Fig [3]. This clearly shows that the spectrum is 6-fold 
degenerate, and more-over, that this degeneracy splits in a fc^. - 3-fold degenerate and a ky- 2-fold degeneracy. This 
is consistent to our analysis, for Ne/N^ = 2/3, N^/Ny = 1/2 which fix q^ =3,qy = 2. Since the GCD{Nx, N^) = 1, 
while the GCD{Ny, N^) = 1, the relative momentum BZ is made oi k^ = and ky = 0,2tt /{Nya) where a is the 
lattice constant in the y direction. In this BZ, the counting of states is 3,4 respectively, while the whole counting 
in the full BZ of the FCI lattice is a 3 x 2 center of mass translation of this in the x,y directions respectively (see 
Figjs]). In the FQH effect, N^/Ng = 1/6 and hence the center of mass degeneracy of g = 6 (= qxqy) and there is no 
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Zero Mode Counting Per Relative Momentum Quantum 
Number of the FQH Haldane Pseudopotential for the 1/3 
Abelian (Laughlin) FQH state 

(Momenta are in units of the smallest allowed 
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FIG. 3; Example of the FQH to FCI mapping for = 2, Nx = 3,./Vy = 4 state. The spectrum of the one-body flat-band 
Chern Hamiltonian ^20^ with nearest neighbor interaction is given on the left. The counting of the low-energy manifold of 
states exhibits clear center of mass degeneracies and the relative momentum BZ of the FCI is a 1 x 2 BZ. The counting per 
momentum sector in this BZ is a folding of the counting of zero modes of the Haldane FQH pseudopotential whose highest 
density state is the Laughlin state. This implies that the universality class of the FCI state is the same as that of the v = !/?> 
state. 



center of mass degeneracy mismatch between the FQH and the FCI. The FQH relative momentum BZ that remains 
after the center of mass degeneracy is removed is a 2 x 2 BZ {GCD{Ns, N^) = 2). The LaughUn 1/3 Abelian FQH 
counting of quasihole states analysis was already presented in Table [l]. In that section, we found that for 2-particles 
in 12 orbitals the quasiholes of the Laughlin state have the following quantum numbers: 2 states in each of the k^, ky 
equal to 0^2n/Ly, 2Tr/Lx,0, and 2^ / L^, 2TT/ Ly momenta, and 1 state in the 0,0 momentum. There is 6-fold center 
of mass degeneracy which multiplies these 7 states to obtain the total of 42 quasihole states. Our theory then implies 
that the counting of the FCI per relative momentum sector is a p^/p = 2- folding of the FQH relative momentum BZ. 
This folding gives a 2 -f 1, 2 -I- 2 counting for the = and ky = 0, 2tt/ (Nyo) momenta respectively, which perfectly 
agrees with the numerical data (see Figjs]). 

For larger number of particles and quasiholes, the cleanest example of the low-energy manifold of states whose 
imprint characterizes the topological order is seen in the (particle) entanglement spectrum of a FCI groundstate. 
The particle entanglement spectrum of the groundstate reveals the physics of the excitations of the system in the 
particular topologically ordered state of the groundstate. It is more robust than the actual energy spectrum because 
the latter can exhibit aliasing effects due to other states at nearby fillings. Since the Hamiltonian used to obtain the 
spectrum is not a "model" Haldane pseudopotential Hamiltonian but a generic Hubbard-like one, adding too many 
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FIG. 4: Example of the FQH to FCI mapping for Ne = 4, = 6,Ny = 6. The one-body flat-band Chern Hamiltonian [20] 
with nearest neighbor interaction for 12 particles in a 6 x 6 lattice has a 3-fold quasidegenerate groundstate from which we 
build the entanglement spectrum shown here. The low-energy manifold (below the dotted line) has 741 states in momentum 
sectors where Nx{ mod 2) = Ny{ mod 2) = and 728 elsewhere. The counting of the low-energy manifold of states exhibits 
clear 3x3 fold (in k^, ky) center of mass degeneracies and the relative momentum BZ of the FCI is a 2 x 2 BZ. The counting 
per momentum sector in this BZ is a folding of the counting of zero modes of the Haldane FQH pseudopotential whose highest 
density state is the Laughlin state. For example, the (0, 0) relative momentum sector of the FCI contains 741 states which 
exactly equals the Px/p ~ Py/p = 2 folding 185 -I- 185 -I- 185 -I- 186 of the FQH Laughlin quasihole counting. Likewise, the (0, 1) 
FCI relative momentum sector contains 728 states, which exactly equals the Px/p ~ Py/p = 2 folding 182 -I- 182 + 182 -I- 182 
of the FQH Laughlin quasihole counting. This implies that the universality class of the FCI state is the same as that of the 
V = 1/3 state. 



quasiholes to the system can, in the energy spectrum, take us to a state at a different filling (with different quasihole 
counting). However, the entanglement spectrum will keep the topological order (filling) of the original groundstate, 
and reveal directly the quasiholes of that state. For example, if we are looking at the energy spectrum of 4 particles in 
36 orbitals on the lattice, it is possible that instead of seeing the spectrum of 4 electrons in a 1/3 Abelian state with 24 
extra-quasiholes, we could, since the interaction is generic Hubbard, obtain the spectrum of 4 electrons in a filling-1/9 
Abelian state or 4 electrons in a filling 1/7 state with 8 extra-quasiholes. However, the particle entanglement spectrum 
of the groundstate with = 12 particles in 36 orbitals (filling 1/3) would reveal the physics of the excitations of an 
Abelian 1 /3 state if we look at the particle entanglement spectrum (described below) of Na — 4 particles versus the 
rest. 

The spectrum of the Chern insulator at 1/3 filling has a quasi-degenerate (3-fold, as our analysis would imply) 
groundstate, at momenta consistent with the folding rule (we wont expand on this rule for the groundstate, as we 
will give the much more complicated quasihole example) . The 3- fold degenerate groundstate of the system contains 
information about the Abelian Fractional 1/3 character of the excitation spectrum. This comes out of the recently 
developed entanglement spectrum '11' T2 which for a single non-degenerated ground state I^P) can be defined through 
the Schmidt decomposition of j^*) in two regions A, B (not necessarily spatial): 
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l*)=E^'^-^'l^^>®l*f> (75) 

i 

where |*/) = (*f |*f ) = The exp(-Ci) and j^-f) are the eigenvalues and eigenstates of the reduced 

density matrix, pA ~ T^t^bP, where p = |^) is the total density matrix. There is no generalization of the Schmidt 



decomposition ( 75 ) to degenerate groundstates but definition of the entanglement spectrum (ES) through the reduced 
density matrix still holds. We want to build a density matrix which commutes with the total translation operators 
-a desirable feature if we want to sort the with respect to the momentum quantum numbers. Such a condition 
is satisfied [TT] by the incoherent summation over the degenerate groundstates p = | X]j Depending on 

whether the A and B regions are real, momentum or particle space, the ES reveals different aspects of the system's 
excitations. If the regions A, B are regions of particles, it can be proved that the particle entanglement spectrum 
obtained by tracing over the positions of a set of B particles gives information about the number of quasiholes of 
the system of A^^ particles and number of orbitals identical to that of the untraced system. In the case of the usual 
FQH, the particle entanglement spectrum of a model state (such as Laughlin, MR, etc) contains an identical number 
of levels as those of the quasiholes, thereby providing a numerical check of our analytic proof [TT]. Away from the 
model states (the Coulomb ground state, for example) the ES may exhibits an entanglement gap [T^JISS]. It separates 
a low-energy structure with perfect quasihole counting and and a high entanglement energy nonuniversal part. But 
a clear and significant gap is not always observed, even for the i/ = 1/3 Coulomb state. 

Surprisingly, for the Fractional Chern insulator, the situation is much better: in a previous work |;8, , we observed a 
clear, large entanglement gap between a manifold of low entanglement energy levels and a manifold of high entangle- 
ment energy levels. Moreover, the counting of the levels below the gap is identical to the counting of quasiholes of Na 
particles in N^Ny orbitals. We show this for a large-size example. Out of the 3-fold quasi-degenerate groundstate of 
12 particles on the Chern insulator on the checkerboard lattice with x Ny = 6 x 6, we build the density matrix and 
construct the entanglement spectrum for Na = 4 particles. The one-body Chern insulator model of [20] presented in 



Eq 78 , with a two-body nearest neighbor Hubbard U — 1 interaction in the flat-band limit is used to generate the 3- 
fold quasi-degenerate groundstates. The entanglement spectrum has a gap (see Fig[4]); the counting of the low-energy 
manifold of states below the dot- ted line is 741 in momentum sectors where Nx{ mod 2) — Ny{ mod 2) = and 728 
elsewhere. The total number below this line (26325) exactly matches the number of (1, 3)-admissiblc configurations 
of 4 particles in 36 orbitals. The counting per relative momentum sector is a folding of that of the zero modes in the 
FQH state, as the complete analysis of Fig[4] shows. 

As said in our previous paper ,8 , we find it very revealing that the FCI groundstates obtained here contain much 
clearer information (large, clear entanglement gap) than the groundstates of the Coulomb interaction in the FQH. The 
entanglement spectrum shows that the groundstates by themselves know of the fractional nature of the excitations 
in the Fractional Chern insulator. The current clean application of the entanglement spectrum also shows that this 
quantity is fundamentally useful towards revealing the physics of strongly-correlated states besides the usual FQH 
model wavefunctions. This is even crucial in this case since, despite several attempts |36l I37j . there is no analytical 
expression of the Laughlin state for the FCI to compare with. 

NUMERICAL PROOF OF EMERGENT SYMMETRIES IN ABELIAN AND NON-ABELIAN 

FRACTIONAL CHERN INSULATORS 

In the current section, we expand our numerical analysis of FCI states to the non-Abelian case and show that 
the counting rule obtained in the current analysis is robust. We perform exhaustive numerical calculations of both 
the energy and the entanglement spectra for several interactions at the specific filling where non-Abelian states are 
expected. The calculations serve as both a numerical check of our emergent translational symmetry principle, as well 
as the first proof of principle that non-Abelian states can be stabilized in a fractionally filled Chern insulator. 

For our one-body model, we pick the Chern insulator on a checkerboard lattice, first introduced in [SKTUlIin]- This 
model already exhibits weak dispersion of the bands (see Figjs]) like several other models [38U40j . But because we 
work in the flat-band limit, this is not essential to our calculation. With A^B being the two sites in the unit cell, the 
one-body, two-band Haniiltonian reads |20) : 

Hi = ^{cl^,clg)hi{k){ckA,CkB)'^ , hi{k) = ^di{k)ai (76) 



24 




FIG. 5; a) One-body model hopping amplitudes in the checkerboard Chern insulator model used in our numerical calculations. 
Direction of the arrow means positive hopping amplitude, b) Dispersion of the model shows that it is an insulator with a nearly 
flat valence band. Our diagonalizations are in the flat-band limit where the valence band is made completely flat, c) Simplest 
Hubbard two-body interaction couples nearest neighbor sites used to obtain the Abelian 1/3 state, d) Simplest three-body 
Hubbard interaction used to obtain the non- Abelian MR Pfaffian state. Equal weight is given to the three configurations e) 
Simplest four-body Hubbard interaction configurations used to obtain the non-Abelian Z3 RR states. Either of the interaction 
1, 2, 3 can be used, all giving a RR groundstate but of different gaps (the largest gap being obtained with the interaction 1). 

where the c?i(A:)'s are: 

dx{k) — Ati cos{(j)) cos{kx/2) cos{ky/2), dy{k) = 4ii sin((/)) sin(fc2;/2) sin(fcj,/2), — 2t2{cos{kx) — cos{ky))+M (77) 

In the original model [20] there is an additional diagonal term -At^ cos (fc^) cos (/c^) - which shrinks the dispersion of the 
bands, thereby making them flatter, but does not matter for the eigenstates; since we are diagonalizing in the flat-band 
limit, we neglect such a term, (j) is the phase factor added to the nearest neighbor hopping, while the parameter M is 
a mass term added in order to drive the transition from a topological Chern insulator (for M = 0) to a trivial atomic 
limit insulator when M — > ±00. The model is always gapped (for ti,t2,4> not vanishing) with the exception of the 
points kx — 0,ky = tt, AI = —At2 and k^ = iTjky = 0, M — 4^2 which are gapless and where the phase transitions 
between the atomic limits AI — > ±00 and the Chern insulator phase occur. For |Af| < 4|t2| the fllled valence band 
has a Hall conductance of 1. The single particle Hamiltonian matrix has the following symmetries: inversion with 
identity inversion matrix hi{—kx, —ky) = hi{kx, ky), as well as (at M = 0) a certain type-of particle hole symmetry 
coupled with a C4 rotation and a mirror operation: azhi{kx, ky)az — —hi{ky, k^). Due to the presence of fractions 
k/2, the model in J20] is not in Bloch form. To render it in Bloch form, we perform the gauge transformation 
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FIG. 6; The spectrum of the three-body Hubbard interaction at filling u = 1/2 shows two quasi-degenerate states at (0,0) 
momentum and 4 states (separated into two quasi-degenerate doublets) at (0,2) (in units of the smallest lattice momentum). 
The procedure and folding from FQH to FCI is also presented. 



CkB Cfcs exp(— ^(fcj, — ky)/2) to obtain: 

hii{k) = 2t2{cos{k^) - cos{ky)) + M (78) 

The inversion symmetry of hi{k) translates into a symmetry of h2{k) given by W {k)h2{k)U (k) = h2{—k) with U{k) 
a diagonal 2x2 unitary matrix with 1^ e~'*^'^"=~'^f on the diagonal. To eliminate the effect of the band curvature, 
and to allow the filled particles to democratically sample the whole BZ, we always work in the flat-band limit of 
a topological insulator. This corresponds to keeping the single particle eigenstates of /i2(fc) but putting the energy 
of the occupied bands to be an arbitrary energy ±i?o where Eq > 0. At the Hamiltonian level, we transform from 
/i2(fc) = E_{k)P_{k) + E+{k)P+{k) to h^^{k) = -EoP_{k) + EoP+{k) where P± are the projectors onto the occupied 
and unoccupied bands. The energy difference between the valence and conduction bands can hence be made large 
without changing the eigenstates of the system. 

We now fractionally fill the flat valence band of this insulator and add interactions. We diagonalize the interaction 
Hamiltonian directly in the filled band, neglecting the conduction band. This is similar to the Lowest-Landau Level 
(LLL) projection in the usual Fractional Quantum Hall Effect. Several interactions we use are: 



-fftwo-body = 2^ nATlB 
<AB> 



(79) 



26 



0.165 
0.16 
0.155 
>, 0.15 
03 0.145 
^ 0.14 
0.135 
0.13 
0.125 



FIG. 7: Low-lying spectrum of the three-body interaction for A^e = 12 electrons and Nx = Ny = 5 

where A, B are nearest neighbor A,B sites. In the current section we also use other many-body interactions such as 
three-body: 

^^throe-body = ^ nAHs^ riBj + '^Sr^-Ai ^^Aa (80) 

<ABi>,<A,B2> 

where Bi, B2 are the two nearest neighbors of the site A and Ai, Ai are the two nearest neighbors of the site B. The 
first of the 3-kinds of four-body interactions used in our exact diagonalization reads (see Fig [5]): 

-ff four-body 1 = ^ ?^Al"-A2"•Bl"S2 (81) 

<Ai,Bi>,<Ai,B2>,<<Ai,A2>>,<<Si,B2>> 

where Bi,B2 are the two nearest neighbors of the sites Ai, and also of the site A2. The other two four-body 
interactions used which we denote i?tour-body 2, ^^four-body 3 can be easily understood from FigjS]. All the numerical 
calculations are performed with t2 — {2 — y/2)/2ti as discussed in [20] • The total translation operators in the x,y 
directions commute with both the single and many-body Hamiltonian and hence the eigenstates are indexed by total 
momentum quantum numbers {Kx,Ky) which are the sum of the momentum quantum numbers of each of the N 
particles modulo {Nx,Ny). The basis states are Jlili 7^ - ...7^ - |0) (we work in the "LLL", and the 7^ -'s are 

the creation operators for a particle of momentum k in the valence band). When acting on the basis states, the 
'^k a ~ ^- a feT- k whcrc u_ ^ ^ IS the a — A, B component of the eigenstate of the occupied band of h2{k) or h2^{k) 
(they have identical eigenstates). Diagonalizing directly in the valence band provides for large numerical efhciency, 
allowing us to go to very large system sizes. The inversion symmetry of the single particle problem is maintained at 
the interacting level, and the interacting spectrum has an exact {Kx,Ky) — >■ {—Kx,—Ky) symmetry which can and 
has been used as checkup. The -fftwo-body has been used to obtain the Abelian, filling 1/3 FCI topological state. 

We now use -ffthrco-body and, borrowing from FQH experience, we look for the MR Pfaffian at filling v = 1/2. We 
diagonalize -ffthroc-body for up to = 14 particles on a A^^; x A^j, = 7x4 lattice, and present the spectrum in Fig[6]. The 
spectrum is separated into a six-fold quasidegenerate groundstate (at momentum (0, 2) we have four states, separated 
into two exactly degenerate doublets) at momenta consistent with our FQH to FCI mapping (in this case we also have 
to use the center of mass degeneracy mismatch factor, which is 2). A clear and detailed presentation of the resolution 
of states in momentum sectors is shown in Fig[6]. Both the number of groundstates and their momenta are consistent 
with the MR Pfaffian state, with a large gap separation to the continuum of states. We now add quasiholes to the 
system, by performing the diagonalization of the 3-body Hamiltonian for A^e = 12 electrons on a A^; x A'j, = 5x5 
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FIG. 8: Entanglement spectrum for the groundstates of the three-body interaction for A^'e = 14 electrons and Nx = 7, Ny = 4 
and Na = 6 particles. The number of states per momentum sector below the dotted line is 8463 in the even ky sectors, 8486 
in the odd ky sectors. This is the counting predicted by our FQH to FCI mapping. 




FIG. 9: Left panel : Low-lying spectrum of the four-body interaction for A^e = 15 electrons and = Ny — 5. The groundstate 
manifold is made of the 10 low energy states in the {kx = 0,ky = 0) sector. Right panel : Entanglement spectrum for the 
groundstates of the four-body interaction for A^e = 15 electrons and Nx = Ny — 5 and Na = 5 particles. The number of states 
per momentum sector below the dotted line is 2050 except in the {k^ — 0,ky — 0) sector that has 2055 states. This is the 
correct counting predicted by our FQH to FCI mapping. 



lattice. This represents the addition of one flux quantum to the system. The spectrum is seen in Fig[7], and we 
observe a low-energy manifold of states (below the dashed line) with 7 states in each momentum sector separated by 
a clear energy gap from a high-energy continuum of states. This is indeed the correct counting for the non-Abelian 
quasiholes of the MR state based on our emergent symmetry principle. We then close the analysis of the MR state 
by looking at the particle entanglement spectrum of the A^g — 14, = 7, Ny = 4 degenerate states when the number 
of particles not traced out of the system is Na — 6. By our theory, this spectrum should show the counting, per 
momentum sector, in a low energy manifold of the quasiholes of the MR state for 6 particles in 28 orbitals . The 
results are presented in Figjs], where we again see that a low-lying manifold of states is separated by a gap form a 
manifold of spurious continuum. The number of states per momentum sector below the dotted line in Figjs] is 8463 
in the even ky sectors, 8486 in the odd ky sectors, the correct numbers of non-Abelian quasiholes predicted by our 
FQH to FCI mapping. 

We conclude that we have presented sufficient proof that the FCI state at half-filling with three-body potential is 
of MR type in the FCI. We now move to the four-body interaction. The Z3 RR state occurs at filling factor 3/5. 
We have performed the calculations up to A^ = 15. Fig. [9] shows the energy spectrum for TV = 15 with a lattice size 
N^ — Ny — 5. We see a tenfold low energy manifold at momenta consistent with our FQH to FCI mapping, clearly 
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separated from the higher energy continuum. The particle entanglement spectrum of this groundstate manifold is 
given in Fig. [9] While all Na sectors display a full entanglement gap, we have chosen the Na = 5 sector. This case 
exhibits an identical counting (2050 states) in all momentum sectors except in the {kx = 0,ky = 0) sector where the 
counting is 2055. Once again, this result is in agreement with our FQH to FCI mapping. 

CONCLUSIONS 

In this manuscript, we have explained the emergent translational symmetry apparent in the low-energy manifold 
of states of a fractionally filled Chern insulator when the insulator exhibits a topologically ordered state. We have 
shown that the one-body GMP algebra obtained at long wavelength in [15 maps the problem of a Chern insulator 
on a Nx X Ny lattice in the problem of a FQH state in the LLL with number of fluxes Ns = N^Ny piercing the torus 
and with a periodic background potential superimposed on the many-body potential already present in the system. 
We revisited the continuum FQH problem and presented a detailed derivation of the many-body relative translational 
symmetries in the system (after the center of mass motion is taken out) presented in a classic paper by Haldane 
[T^ . We then showed that the BZ of the FCI is just a folding of the one for the FQH, and pointed out the direct 
connection between the center of mass degeneracies in the FCI and those in the FQH. This FQH-FCI map allows us 
to compute the number of states that arise per momentum sector of the FCI if the FCI is in the same topological 
phase as the FQH. We showed how to obtain, in the FQH effect, the counting of quasihole states per 2D relative 
momentum sector for a large series of states exhibiting a generalized Pauli principle [H]. These states include the 
Laughlin, MR and RR series, as well as the Jack polynomials states, the Haffnian|41|. and many others, including 
extensions to spin-unpolarized states [^l43j . The counting construction makes use of only combinatoric methods and 
involves counting partitions. We showed how to use the combinatoric counting rule and the FQH-FCI map to obtain 
the number of states per momentum sector in both the energy and the entanglement spectrum of the FCI at filling 
1/3. We then proceeded on a comprehensive numerical analysis of the spectra of Chern insulators with three and 
four-body interactions. We provided a proof of principle that Moore-Read states and Read-Rcazayi Z3 non-Abelian 
states can be stabilized in Chern insulators and provided a numerical check of our FQH-FCI mapping. 
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APPENDICES 

Appendix 1: Details of the Girvin-MacDonald-Platzman Algebra 

In this appendix, we prove in detail that 

E,,ni* "2 "2* "12 _ "1* "2 "2* m2 Sri l('n-7/;- - 7/IY7-'|F"l'"2 (oo] 

^a,k a.,k+q^l3.,k+q^l3,k+w+q ",3,fe+tu "a.fc+tu "a,fc+«;+g ~ 2 ^^^O'^ij 

"2 

The expansion (double index means summation): 

<,k+q = K.k + <ltdtK,k + \<li<l]didjul^, (83) 



' 2 ' 

^ni ,m2 



gives: 

Km + ^i^»<'fc + h'^iWjdtdju2^){u'^l + [q, + Wi)diU^l + \{qi + w^){qj + = 

= iSm,n2 + qiUl%diul\ + h^qjKSdrdju'^^) ■ {Sn2.m2 + + W^ufJ^diuJl 

+ \{q^ + w,){q, + w,)u-%d,d,u^l + q^uj^d^u-^ + g,(g, + w,mu-%){d,u^l) + \q.q,u^:Ad,^7J:) " 

+ h^<ljK]kdidju'^X -iq^w,a^ 13) 
The zeroth order term in q cancels. We now group the first and second order terms. The first order terms are: 

[qi + w^)u^p]ldiuj^f^ + q^^^{uJ].u^]^) -{q^w,ai^ P) = 

{qt + 'Wi)u^'^^diujl - (q O w, a ^ /3) = 

(summation over a, j3 was implied. Hence the first order term is zero. 
The second order term is (remember double index summation): 

\{q^ + w,){q, + W,)uf*d,d,uJX + qM, + wm^7;k)id,uJX) + 

H(i^<io'^Tk^^^j^}S + 9^(9* + ^i)'^l%^o^l]khls^kdiuJl + q]qtKSi9jul''^k>^,l^iu''^S + 5*9j<!fc5»5j<l- - 

— (gr w, a f^- /3) = 

+q,q^<S{^,u2,)u^l^,u-^% - w,w.u;^:{d,u;^,)u::xd.ui% + ^{q.q, - w.w,)u-jp,d,uZ = 

= {q^Wj - w^qj){{diu'p]^){djuj^^) - ul]l{djul\)up'ldiU^l) + 
+ i(9.q, - w,w,){u^p,d,uf; + ul]ld.,d,u2, + 2{d,uf;){d,u^^)) + 
+<lm{<%^,ul^,k>^l^^^-% + u2l{d,ul^,)u-%d,u^X) - ^^^A^;S(d,u;^k>7kd,ul% + uf;,(d,u}^^,)u2ld,u 
We now observe several simplifications. The last line of the above equation vanishes due to the identity: 

The second but last line of Eq[86] also cancels because of the identity: 

when Aij is any symmetric matrix. We are left with 



n2 "'a,k"'a,k+q"-l3,k+q"-l3,k+w+q P ,k ^ 13 ,k+w^a,k+w"'cx,k+w+q ~ 

{qiWj - w^qj){{^^uf *){^juJl) ~ ul\l{djul^,;)uf^ld^uj^,;) = 
= \{q,w, - w,q,mu;]:){d,uj^,) - ul%{d,u2,)u}%d,uJl - (z O 3)) 
By denoting the non-Abelian Berry potential 

we find that 

where F^/^""' = 9,^;^^"^ - d^A^''''^ - 
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Appendix 2: Relation Between the Density Algebras and Chern Number 



The projected densities in the commutator algebra in Eq 60 are required to not commute if the topologi- 
cal insulator is to have a non-zero Chern number. Taking the trace of the commutator in Eq 60 with the 
P-q-w = Y.k ^"ak-q^w'^'aklT^ 1^) (^1 Ik-q^w obtaiu, to sccoud ordcr in q,w (in the long- wavelength limit): 

[p^^p^]p_^_^^^'-{q.,u,,~w,q,)-^ ^ i^i'""^(fc)7r^|0)(0|7r (92) 
J-ak-q-w'^ak = ^mi.n2 to zcroth Order in q, w: 
[Pq.s.,P^JP-q.^-n.,y = lq.Wy-^ ^ i^r,""nfc)7r^ |0) (0|7r (93) 



^3/2 



We now take the trace of the above to obtain 

2 x/]vr 



. Tr{[p,^,,p^^y]p^,^,^^^y) = ^^F™^™(fc) = C (94) 

^ k,m 

with Qx = ^ , Qy = the lowest value for the momenta one can obtain on the lattice. Hence the Chern number can 
be expressed as the long- wavelength limit of the commutator of the projected density operators in the lower sub-band. 

Appendix 3: Relation Between the Position Operators and the Chern Number 

An alternative way of presenting this result is by using the position operator in the X direction; 

where by Akj we mean Afc • j. The projected density operator in the lowest band is: 

P^P=^T. OT+Aqllk |0) (0| iT+Aq = PAq (96) 



Hence the projected density operator for Aq — 2t: /N,^. equals the projected position operator. There is an alternate 
expression for the projected position and density operators. They are translational operators in fc-space: 

PXP \k, n) = PA, \k,n)^ K'-Aq.A Ifc - "i) « e-^-^-^^"' |fc ~ Aq, n,) (97) 

For one occupied band, the density or position operators just translate the band at different momentum, whereas for 
more than one occupied bands, the translation is accompanied by a rotation. We now see that: 

2 n'^^'^ 1 

-^J^Tr{[PXP,PYP]X-^Y-^) ^ ^ 5]i^r,"(fc) - C (98) 

k,m 

Appendix 4: Chern Number of a Landau Level 

We outline how a simple calculation for the Chern number of lowest Landau Level. The wavefunctions for the 
landau levels in the continuum do not have two momenta kx,ky, and hence the momentum expression for the Chern 
number cannot be used. We obtain the Chern number as: 

C = 4 Ajk--.o ^n[PXP,PYP]X-^Y-^] (99) 
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where in the LLL we have the expressions 



PXP = ^Pe*^»^P = Pei^-(^+^)/^)p, PYP = -^Pe'^^yp = Pe^''(^-^)/2)p (100) 



LLL wavefunctions tp{z) have the property that Il'ip{z) = where 11 = IIj; + iILy, IL^, Ily being the canonical momenta 
in the presence of a magnetic field. We then have: 

p^za+zb p ^ pgi^(na-ntb)gi^(J<-t(,_ifa) p {101) 

where we have used z = {il'^/h){Il — K), K = + iKy, and the fact that 11 commutes with K. We now need 
to disentangle the exponent containing n,nt to obtain: e'^Cna-ntb) ^ ^-i'^uH ^i'^na ^abi'' / h\ gjj^^^g e'^nop = 
Pe-^Tn+b = p^ we have 

p^za+zbp _ ^abF /h'^ p^i^{K^b-Ka) p (102) 

. This allows us the computation of the terms involved in the Chcrn number: 



[PXP, PYP] = le-^^^^[Pe'T^=^-ff!/p,Pe-4^>/-ff^p] (103) 
Moreover, again after projection to the LLL, and after some algebra, we find: 

Tr[[PXP,PFP]X-iy-i] = ^e-s^^"^ x 

J'j.^^-j^^^g-i^A^Ky p^i^A^Ky p^-i^AyK^ p^i^AyK^ _ pj^^^ ^i'^AyK^ p^-if^AyK^ p^i'^A^Ky pg-i4'^-fl<|4) 

^2 . j2 

We now observe that, if we write down P = J cPr liplr)) {ip{r)\ then the quantity e^'^'f^^'^^y Pe^'fi^^^^ is 

^2 

again a projection operator P due to the fact that e~*"S"^"=^'' are translation operators. We find that: 
Trfe-'T^^-^yPe'T^-^j/Pe-'x^s/^^Pe'T^i/^-] = V. We hence obtain: 

^ h'^ 1 i2 A^+Ag t2 A^A„ ,2 A^ 

C = --5- lim e~^^~(e~^^~ - e^^~) = 1 (105) 

P A^,Ay^O A^Aj, ^ ^ ^ ' 



Appendix 5: Center of Meiss Degeneracy Mismatch 

The center of mass degeneracy in the FQH is g = Ns/GCD{Ns, Ne) (with Ns = N^Ny) while in the FCI is 
q^Qy = {N^/GCD{N^,N^)){Ny/GCD{Ny,Ne)). Their ratio is 



GCD{N^,Ne) ■ GCD{Ny,N^) 



(106) 



q^Qy GCD{N^Ny,N^) 

We can easily prove that this is an integer by applying the decomposition of a number in primes Pi'. = YiiPt*' 
Ny = Y\iPi \ Nf, = Hi Pi' where ai,l3i,9i arc all integer powers. The ratio then becomes: 

q _ TT min(ai,6li)+min(ft,ei)-mm(ai+/9i,6li) (107) 

q^qy ~ i ' 

We can now analyze the different cases to prove that all powers of pi are positive in the above expression. Since 
our analysis is valid for any i, we drop the index. We have the cases l)a<0,/3<9,a + l3>6m which case 
min(a,6') + min(/3, 6*) - min(Q; + 13,6) = a + /3 ~ 6 > ih 2) a < 0, (3 < O.a + B < 6 in which case mm{a,0) + 
min(/3, 9) - mm{a + /3,0) ^ 0; 3) a > 9, P < 0, a + (3 > in which case min(Q!, 9) + min(/3, 9) - mm{a + 13,0) = (3; 4) 
Oi>9,(3 < 0,a+/3 < 6'inwhichcasemin(a,6')+min(/3,e)-min(Q;+/3, 6') = 0~a >0;5)a>9,P> 9,a+f3 > eiin which 
case min(Q:, 0) + min(/3, 0) — min(Q: + (3,0) = 0. We see that in all cases possible, min(a, 6) + min(/3, 9) — min(a + /3, 9) 
is a positive integer which makes q divisible by qxqy 
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